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Abstract. In this paper we develop an existence theory for small amplitude, steady, two- 
dimensional water waves in the presence of wind in the air above. The presence of the wind is modeled 
by a Kelvin-Helmholtz type discontinuity across the air-water interface, and a corresponding jump 
in the circulation of the fluids there. We consider both fluids to be inviscid, with the water region 
being irrotational and of finite depth. The air region is considered with constant vorticity in the case 
of infinite depth and with a general vorticity profile in the case of a finite, lidded atmosphere. 
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1. Introduction. Understanding the precise means by which the wind creates 
surface waves in the ocean is both a central problem in geophysical fluid dynamics and 
a famously difficult one. From common experience, of course, the fact that wind blow- 
ing over quiescent water will lead to the generation of waves is fairly obvious. Indeed, 
the basic physical process appears straightforward: the jump in tangential velocity 
across the air- water interface causes instability, creating growing modes, which even- 
tually stabilize to become traveling waves. Were this the case, we might hope that 
the salient features of the system might be captured by the Kelvin-Helmholtz (K-H) 
model (cf., e.g., [!)]). Intriguingly, this appears to be untrue. Including the effects of 
surface tension, the K-H model predicts that the speed of the wind must be above 
660 cm/sec in order to excite waves, which is an order of magnitude larger than 
what observation suggests (cf. [15]). The discrepancy indicates that there are crucial 
components of the system that has been overlooked by the K-H viewpoint. 

The search for these missing features has led to a number of competing models 
for the time-dependent wind-driven generation of water waves with perhaps the most 
influential being that of Miles (cf. [17, 18]), which is based on the existence of a 
critical layer in the continuous wind shear profile above the surface wave (we will 
discuss some aspects of this model in §1.2 below). 

The authors will address the time-dependent problem in a later work, but in the 
present paper we begin by consider the steady problem. That is, we investigate the 
question of existence of traveling waves in a two phase air-water system. We endeavor 
to do this in such a way that we can view the waves as the eventual byproduct of 
wind blowing over water, though we shall remain agnostic as to how exactly that 
generation took place. For us, this means that the circulation as a function of the 
streamlines should be discontinuous over the air-water interface. Since the circulation 
is conserved by the flow, this is a necessary condition for the traveling wave to be 
dynamically accessible from an initial configuration with laminar flow in the air and 
water with a jump in the tangential velocity over the boundary. 

In this section, we shall first describe the basic framework, then make an informal 
statement of the results. The precise theorems will be given later in the text. 



tCourant Institute of Mathematical Sciences, New York University, New York, N.Y. 10012 
twalsh@cims.nyu.edu 
§ obuhler@cims.nyu.edu 
^shatah@cims. nyu.edu 



2 S. WALSH, O. BUHLER, J. SHATAH 





Y 

Y = £ 




fii(t) 

Y = V (X,t) 




' X 



Fig. 1.1. The fluid domain U(t). The wave is assumed to propagate to the left with speed c. 

1.1. Eulerian formulation. Consider a two-dimensional surface wave in the 
ocean, with an accompanying wave in the atmosphere above. Fix a Cartesian coor- 
dinate system (X, Y) so that the X-axis points in the direction of wave propagation, 
and the Y"-axis is vertical. We assume that the bed of the ocean is flat and occurs at 
Y = —d, while the interface between the wave and the atmosphere is a free surface 
given by the graph of a function r\ = r](X, t). We are particularly interested in the case 
where r\ is periodic in X, for fixed t. While these are not physical (they have infinite 
energy), they are the type typically studied in connection with linear stability. Not 
coincidentally, they are also more amenable to mathematical analysis since there is a 
gain of compactness in the X-direction. We may normalize rj by choosing the axes so 
that the free surface is oscillating around the line Y = 0. The atmospheric domain can 
be thought of as either unbounded or bounded in Y. The unbounded regime models 
the situation in which the characteristic horizontal length scale is vastly smaller than 
the vertical length scale. On the other hand, if the dynamics of the wave away from 
the interface are deemed of secondary importance, a common practice is to take the 
air region to be bounded above by a rigid lid at Y — I. All told, the fluid domain at 
a given time t is 

fi(i) = fii(t) Ufi 2 (*), 

where 

ni(i) := {(X, Y) e K 2 : V (X, t) < Y < 1} 

is the air region (£ = +oo in the unbounded case) , and 

fi 2 (t) := {(X,Y) G K 2 : -d < Y < rj(X,t)} 

is the water region. The total width of the channel is thus W := £+d. In what follows, 
we shall assume that t is fixed at the outset while d is to be determined during the 
solution procedure. We point out also that we are not including the air-sea interface 
{Y = rj(X,t)} in the fluid domain, and thus Q(t) is not connected. 

Let u = u(X,Y,t) and v = v(X,Y,t) denote the horizontal and vertical fluid 
velocities, respectively. Let g = g(X, Y,t) > be the density. We assume that 
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velocity field is incompressible 

u x +v Y = 0, in 17(f). (1.1) 

Conservation of mass is enforced by stipulating that the density of each fluid particle is 
likewise preserved. For an inviscid fluid, this is equivalent to stating that the material 
derivative of g is zero: 

Qt + ugx + vqy = 0, in 17(f). (1.2) 

In this paper, we exclusively consider the case where the density in the air and water 
regions is constant, that is, 

g(t,X,Y) — PairXnU) (t) {X, Y) + PwatcrXf2(2) (t) (X, Y) , 

where Xn(^(t) 1S the characteristic function for the fluid region 17W(f), i = 1,2, and 
Pair, Pwater are the given densities of water and air, respectively. Thus (1.2) will always 
be satisfied. Note that because 17(f) does not include the air-sea interface, the above 
equation will hold in the classical sense. 

The momentum equations for non-diffusive heterogeneous fluids are the Eulcr 
equations, 

eut + euux + evuy = -p x , mn{t (13) 

gv t + guv x + QVVy = -Py — gg, 

where P = P(X, Y, f) is the pressure and g is the gravitational constant. Again, given 
our choice of g, one may alternatively view this as being satisfied in fl^(t) and fl^ 2 \t) 
separately, with g taking the appropriate constant value. 

Let I[t) := dQ,\(t)C\dQ.<2(£) denote the interface between the air and water regions. 
For simplicity, we shall use the convention that the restriction of any quantity defined 
on Q(t) or 17(f) \ I(t ) to the subset 17j(f) is denoted by a superscript (i). Thus, for 
example :— u|n ( (t)' Similarly, we define the jump operator 

H := (■)%-(■)%■ 

The dynamic boundary condition states that, ignoring the effects of surface ten- 
sion, the pressure must be continuous across the interface. Stated in terms of the 
operator [•], this is simply the statement that 

[PJ=0, on Y = r)(X,t). (1.4) 

To couple the evolution of the boundary to that of the flow, we impose a kinematic 
condition. More precisely, we suppose that 1(f) is a material interface. This is enforced 
by requiring that the normal velocity of the boundary agrees with the normal velocity 
of the fluid. Since we are assuming a graph geometry for the free surface, we can 
express this quite explicitly: 

v = r]t+ur)x, on Y = r)(X,t). (1-5) 

Similarly, in the lidded regime, both upper and lower boundaries are unmoving, we 
must have 



v = 0, on Y = —d. 



(1.6) 
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v = 0, on Y = L (1.7) 

When is unbounded, the bed condition (1.6) remains valid, but (1.7) must be 
replaced with 

(u, v) — > (uooj 0) as Y — > oo, uniformly in X and t, for some Moo G K.- (1-8) 

Traveling wave solutions to (1.1)-(1.6) are those where, for some wave speed c > 0, 
the change of variables 

x = X — ct, y = Y, 

eliminates time dependence. Physically, this means that viewed from a frame of ref- 
erence moving with fixed speed c in the direction of propagation, (u, v,g, rj, P) all 
appear steady. Reusing notation, from here on we shall simply identify (u, v, g, rj, P) 
with their stationary profiles. Observe that periodicity of the traveling wave is equiv- 
alent to periodicity of the steady quantities in the x-variable. We shall therefore 
require that (u, v, g, rj, P) are L-periodic in x, for some L > 0. 
In the moving frame (1.1)-(1.3) become 

U X +Vy = 0, 

g(u — c)u x + gvuy = —P x , in fi, (1-9) 

g(u - c)v x + gvvy = -P y - gg, 

where Q is the steady domain. The kinematic and dynamic conditions for the lidded 
problem are likewise 

= 0, ony = l, 

- ?' , on y = -d 
= (u-c)r) x , on y = r](x), 
= 0, on y = r)(x). 

The unbounded atmosphere case differs only in the condition at y = oo, where we 
require 

(u, v) — ¥ (uoo,0) as y — > oo, uniformly in x for some Uqo € M. (1-H) 
Recall, also, that we have chosen our axes so that 77 oscillates around the line y = 0: 




L/2 

r](x)dx = 0. (1.12) 

-L/2 



We note in passing that the steady Euler equations have an important consequence 
for the Eulerian-mean momentum flux defined as 

L/2 

F E (y) = \ J p(u-c)vdx. (1.13) 

-L/2 

Physically, the function Fe(v) is the mean upward flux of horizontal momentum across 
a line of constant altitude y. The steady Euler equations imply that dF E /dy = and 
hence F E does not depend on altitude y. This constant-momentum-flux result is also 
obvious from a physical point of view: any vertical variation of F E would imply a 
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time-dependent accumulation of horizontal momentum within some altitude range, 
which would be incompatible with the assumption of a steady flow field. Moreover, 
the boundary condition v = at the rigid lower boundary actually implies that Fe = 
there and therefore Fe = throughout the domain. This condition on Fe must be 
satisfied by all solutions to the steady equations. 

Finally, another physical quantity of interest in the context of wind-driven waves 
is the mean horizontal drag force exerted by the air on the water across the interface 
y = r){x). This drag force equals minus the Lagrangian-mean momentum flux 1 across 
the undulating air-water interface, which is 

L/2 

drag force = — Fl = j- J rj x (x)P(x, r)(x)) dx. (1-14) 

-L/2 

It is easy to show by integrating the steady Euler equations over a control volume 
with lower boundary y — T]{x) and an upper boundary of any constant altitude y 
above the interface that 

Fe = Fl = (1.15) 

holds for steady waves, i.e., the Eulerian momentum flux equals minus the drag force 
across the interface, and both are zero for steady waves. Conversely, a nonzero drag 
force is incompatible with a steady flow, and this observation lies at the heart of 
Miles 's theory for wind-driven water waves, which is briefly discussed in the next 
section. 

1.2. Critical layers and Miles's theory. From (1.9) and (1.10) it is clear 
that points where u = c are of special importance. When this occurs, the relative 
horizontal velocity appears to vanish in the moving frame, meaning that the horizontal 
velocity of the fluid particle at that point matches that of the wave. This scenario 
we refer to as stagnation. Note that this differs from standard usage, since for us 
only the relative horizontal velocity needs to be zero. In the classical theory of steady 
water waves, stagnation is closely associated with a loss of regularity stemming from 
the degeneracy of the governing equations. The most well-known instance of this 
phenomenon is the so-called extreme waves of Stokes (cf., e.g., [2]). Stagnation points 
play an even more central role in the time-dependent theory of wind-driven water 
waves. Note that when a flow is laminar, i.e. it is of the form (u,v) — (U(y),0), 
the critical points will arrange themselves in horizontal lines y — y c , say, such that 
U(y c ) = c. These lines are called critical layers. The central thesis of Miles's theory 
is that the presence of a critical layer in the air region provides a mechanism for the 
wind to impart horizontal momentum on the water via a nonzero drag force in (1.14), 
and it is precisely this drag force that is responsible for creating linear instability at 
slower wind speeds than predicted by K-H. (cf. [17]). 

Specifically, if the vertical momentum flux Fe at the upper domain boundary is 
zero, which is consistent with a lidded domain, and if there is only a single critical 
layer, then linear wave theory for weakly unstable waves predicts that Fe has a jump 
discontinuity across y = y c such that the drag force on the water is proportional to 
the ratio (—U"/U') evaluated at the critical level y = y c (cf. [17, 9]). This makes 



1 Eulerian and Lagrangian momentum flux definitions and their role in wave dynamics are dis- 
cussed in detail in [5] 
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obvious the crucial role played by a nonzero value of the second derivative U" {y c ) at 
the critical layer, without which there could be no wave growth induced by the critical 
layer. Indeed, if U"(y c ) is assumed to vanish, the critical layer is neutered or inactive: 
the momentum flux Fe is then continuous across y = y c and the critical layer plays 
no part in the generation of waves. 

Now, for the purpose of studying steady waves one must either rule out criti- 
cal layers a priori, or one must allow only inactive critical layers by requiring that 
U"(y c ) = 0. Indeed, the existence theory we present in §6 for the constant vorticity 
unbounded atmosphere case allows for these inactive critical layers. 

The majority of the present paper is concerned with the lidded atmosphere case, 
and in that we regime we restrict our attention to waves without stagnation and 
hence without critical layers. The alternative approach of studying the existence of 
small amplitude traveling water waves with stagnation points in various regimes has 
recently been considered in several works (cf. [21, 10, 11, 16]). In each of these papers, 
the authors prove that there exists a (local) curve of non-laminar flows bifurcating 
from a background laminar flow with a critical layer. Though it is not always stated 
explicitly, in every case this is done under the assumption that the critical layer occurs 
only at an inflection point of U and therefore U"(y c ) = 0. If a similar restriction is 
taken, we believe it would be possible to generalize the results (in lidded domains) 
given here to allow for critical layers in the background flow. As the above discussion 
makes clear, however, in the context of wind-driven waves these types of profiles are 
not considerably more interesting than those without stagnation. 

1.3. Stream function and circulation. Assuming an absence of stagnation, 
(1.9) ensures that we may define a function ip — ip(x,y) by 

ipx = -y/ffv, ijjy = y/g(u - c), in fl. (1-16) 

This tp is known as the (relative) pseudo-stream function for the flow. Here we have 
the addition of a factor of g to the typical definition of the stream function for an 
incompressible fluid. This is meant to account for the inertial effects of the density 
variation (cf. [22, 23]). Observe also that the restriction u < c throughout the fluid 
becomes the requirement: 

ip y < 0. (1.17) 

It is evident from (1.16) that ip is indeed a (relative) stream function in the sense that 
V J -?/' is collinear with the vector field (u — c,v). In other words, the level sets of ip, 
which we call streamlines, coincide with the standard definition of streamlines in the 
literature. 

The kinematic condition in (1.10) implies precisely that the free surface, bed, and 
lid are each level sets of ip. As (1.16) only determines ip up to a constant, we may 
take ip = on the upper lid, so that ip = —po on y = —d, where po is the (relative) 
pseudo-volumetric mass flux: 

rt 

Pa ■= V Q(x, y) [u(x, y) - c] dy. (1.18) 

J-d 

It is straightforward to show that p is a (strictly negative) constant, i.e., it does not 
depend on x (cf. [ ]). Physically, po describes the rate of fluid moving through any 
vertical line in the fluid domain (and with respect to the transformed vector field 
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s/q{u — c, v).) Notice that, although wc shall allow u, v, and p to have discontinuities 
across X, we assume that the pseudostream function is of class C(ft). This is not an 
assumption: because ip is defined by (1.16) only up to a constant in each Q( l \ we may 
without loss of generality take it to be continuous across the interface. 

The conservation of mass in (1.9) implies that Vg is orthogonal to the velocity 
field, and hence we may let p : [—po, 0] — > R + be given such that 

Q(x,y) = p(-il>(x,y)) (1.19) 

throughout the fluid. We shall refer to p as the streamline density function, though 
one may alternatively view it as the Lagrangian density. In the ocean, one typically 
has that density is increasing with depth, meaning that the lighter fluid elements 
rest on the heavier ones. Indeed, several physical mechanisms work independently 
to enforce this situation: gravity naturally leads to increased salinity near the bed, 
while temperature decreases substantially as one moves deeper into the ocean where 
the effects of the sun's heating are attenuated. However, even near the surface, water 
is on the order of 1000 times as dense as the air. The variations in the density are 
within the air and water region are nowhere near as great as that across the interface, 
and so we shall suppose that g is constant in each region, 

Q ^ — Pair; Q ^ — Pwatcr- 

Stable stratification in this case simply means p^ < p wa ± eT . 

Conservation of energy can be expressed via Bernoulli's theorem, which states 
that the quantity 

E:=P+^((u-c) 2 + v 2 )+gpy, 

is constant along streamlines (see [22] for an elementary proof.) If we evaluate the 
jump of E on the interface, we may use the dynamic boundary condition to express 
the pressure in terms of (u, v, n), which gives rise to the following 

[|W| 2 ] + 2.g M (77 + d)=Q, on y = r)(x) (1.20) 

where the constant Q :— 2 \E + gpd\ gives roughly the jump in the energy density 
across the free surface of the fluid. We treat Q as our parameter of bifurcation. 
By taking the curl of the steady Euler's equations, one arrives at the identity 

{-Aip, lu} = in n, 

where {•, •} is the Poisson bracket, and w is the scalar vorticity uj :— v x — u y . Under 
the assumption that u < c throughout the fluid, this allows us to conclude that there 
exists a single- valued function 7, called the vorticity strength function, such that 

u){x,y) = j(ip(x, y)) in fl 

The final ingredient in our model is a condition on the circulation on each stream- 
line, namely that the circulation in the air and water regions do not coincide. This 
is meant to ensure that the waves we construct are dynamically accessible from an 
initial configuration of a shear profile wind blowing over water. Since the circulation 
will be constant along a streamline, we define the function r = T(p) defined by 



T{p):= [ (u(x,y),v(x,y))-dx = r(p), for Po < p < 0, (1.21) 
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here S^(p) := {(x,y) € ft : —L/2 < x < L/2, ip(x,y) — — p}. For us, then, the pres- 
ence of wind means that there is a discontinuity in T at the streamline corresponding 
to the air-water interface. 

1.4. Informal statement of results. Now that the basic objects have been in- 
troduced, we can give a brief summary of our results. The precise theorem statements 
are presented later. 

(i) Consider the existence of steady wind-driven water waves with a lidded 
atmosphere and irrotational flow in both the air and water and without stagnation. 
Fix the period L, density jump [p], (pseudo) volumetric mass flux in the water region 
Pi, the (pseudo) relative circulation on the lid r rc i, and the height of the lid t. Then, 
if a certain compatibility condition is met (3.2), there exists a family of laminar flows 
that are exact solutions. Moreover, if a size condition is satisfied (ILBC), there is a 
curve of small amplitude (classical) exact solutions bifurcating from this family. See 
§3 and Theorem 3.1. 

(ii) Consider the situation as in (i), but where the flow in the air region is 
rotational. There is a corresponding compatibility condition relating 7, r rc i, p\ and 
I (4.2). If it is satisfied, there exists of a one-parameter family of laminar flows, each 
with (pseudo) relative circulation r rc i. Moreover, under a certain local bifurcation 
condition (LBC) (or size condition (4.9)), there is a curve of small amplitude (classical) 
exact solutions bifurcating from this family. See §4 and Theorem 4.1. 

(iii) In the unbounded atmosphere regime, we fix the depth of the ocean d and 
consider the existence of waves where the water region is irrotational and the wind has 
constant vorticity 70. There is a family of laminar flows, parameterized (essentially) 
by the circulation at y = +00, and, under analogous bifurcation conditions ((5.2) for 
70 = and (6.1), otherwise), there are curves of small-amplitude (classical) exact 
solutions bifurcating from this family. See §5-6 and Theorem 5.1 and Theorem 6.1. 
Note that these results do not require an absence of stagnation. In fact, when 70 < 0, 
there will be a critical layer in each of the background laminar flows. 

Notice that in each of (i) and (ii) there are hypotheses relating r rc i, I, 7, and p\. This 
is to be expected, in fact they can be viewed as a consequence of Stokes' theorem. 
In our work, we elect to fix po and pi, as well as t in the lidded case. When the 
atmosphere is irrotational, these choices determine r re i by (3.2); if the atmosphere is 
rotational, we take 7 to be fixed, and define r re i according to (4.2). These choices are 
of course arbitrary and one can instead choose r re i and use (4.2) to define 7 and £. 

Let us now briefly discuss the place of these results in the existing literature. The 
bifurcation theory techniques that we employ have a long history in the study of steady 
water waves. For the lidded regimes (points (i) and (ii)), we consider a reformulation 
of the problem in semi-Lagrangian coordinates, which has been used in a number of 
works, notably Amick-Turner (cf., e.g., [19, 20, 1]) and Constantin-Strauss (cf. [7]). 
Our approach follows the latter in relying on elliptic estimates rather than variational 
techniques. More precisely, the method we employ can be viewed as an adaptation 
of that in [22] to the case of a non-continuously stratified fluid in a channel, and 
with additional considerations involving the (pseudo) relative circulation. A similar 
problem was considered by Amick and Turner in [3], ignoring the important issue of 
the circulation. They, however, were primarily interested in the solitary wave case, and 
so developed the periodic existence theory only in order to obtain solitary waves as a 
limit as the period goes to +00. As a consequence, this requires them to make certain 
assumptions on 7 (otherwise the limiting wave will not be irrotational and quiescent 
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at x = ±00); we do not impose any such restriction. It should be pointed out that 
the results of Amick-Turner, as well as those of Constantin-Strauss, are global, in 
the sense that the bifurcating curve of solutions is extended to include waves with 
finite amplitude. We believe that similar result is possible in this case, since the basic 
ingredients (mainly good Schauder-type a priori estimates for the elliptic equations 
involved) are available. This is an interesting and important question, but beyond 
the scope of the preliminary investigations here. 

Previous work on the infinite atmosphere case is comparatively sparse. In the 
applied literature, this is simply because the lidded regime is seen as an adequate 
idealization of the infinite atmosphere: so long as the wind curvature is evanescent 
and the air density is constant there is no mechanism for propagation of waves to 
or from vertical infinity, and hence the dynamics far away from the air-sea interface 
decay exponentially with altitude and are not thought to be particularly relevant. 
Mathematically, of course, removing the lid results in a loss of compactness, which 
introduces some potentially serious difficulties. Nonetheless, we include as a simple 
application of our machinery a mathematical treatment of the infinite atmosphere 
regime in the case where the vorticity is constant in the air region. 

2. Formulation. 

2.1. Stream function formulation. The relevance of the pseudo-stream func- 
tion and the vorticity strength function to the existence theory stems from the identity 

-A-0 = w, 

which, recalling the definition of the vorticity strength function leads to the semi-linear 
elliptic equation 

- A0 = 7(0), in n. (2.1) 

One important consequence of the lack of stagnation, or, more accurately, the absence 
of active critical layers, is that the Euler system can be recast as the above scalar 
problem. 

Next, we note that the (steady) kinematic condition in (1.9) guarantees that the 
interface I is a streamline. That is, 

I = {0 = -Pi}, 

for some p\ > po- The difference between p\ and the value of tp at the top boundary 
of the domain gives the (pseudo) volumetric mass flux in the air region. In the lidded 
case we thus take pi to be some fixed value and let ip\ y= e — 0. When the atmospheric 
region is unbounded, however, we let ip\x = 0, and thus ip — ¥ +00 as y — > +00. 

Taken together, the considerations of the preceding paragraphs imply that ob- 
taining solutions to (1.3)— (1.21) with a lidded atmosphere for a given p and 7, is 
equivalent to solving the following problem: Find (ip, T], Q) such that ipft < in fiW, 
and 

A0 + 7(V>) = 0, in ft, 

[|VVf] + 2g [p] (r, + d) - Q = 0, on y = r)(x), 
< -0 = on y = £, (2.2) 

ip = -Pi, on y = r)(x), 

, -0 = -Po, ony= -d. 
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For any such solution, the relative circulation will then be given by (2.3). We shall 
forestall a detailed discussion of the unbounded atmosphere case until §5 and §6 

Note that we have defined SI so that it does not include the interface. Thus 
(2.2) can be thought of as two separate elliptic problems in the domains SI 1 - 1 ' and 
Sl^ 2 ), that must then be matched along the interface according to the jump condition 
(1.20). The advantage of this view point is that, while this a free boundary problem, 
the coefficients of the elliptic equation are all smooth. Another way to proceed is 
to pose (2.2) in a weak sense on O, with the jump condition being represented as a 
measure supported on I. This conceptualization allows us to understand the matching 
procedure in the more conventional framework of elliptic problems with non-smooth 
coefficients, for which there is a great deal of theory. As stated, we will be taking the 
first view — namely, that (2.2) is two elliptic problems matched at the interface — but 
occasionally will make use of the second view to derive some compactness properties 
of the corresponding operator. 

Finally, let us discuss the circulation for the reformulated problem. Since S^(p) 
is a level set of ip, and we have 

(«,«) = -^V^ + fc.O), 

the circulation in the air is given by 

r(p) = Lc- f —L^iv^du 1 , 
J.y(p) vp(p) 

where H 1 denotes one-dimensional Hausdorff measure, which is equal to the arc length 
of the interface. It will be more convenient to consider the quantity 

r re i(p):=i/ \V^\dU\ (2.3) 

which we call the (pseudo) relative circulation. T and r rc i arc then related according 
to the equation 

T rcl (p) = ^p(p)(c- r ^). (2.4) 

The advantage of considering r rc i in place of T is merely that it is simpler to express 
in terms of ip, while being equivalent for a specified p, c, and L. 

2.2. Height equation formulation. The main difficulty that remains in (2.2) 
is that the domain SI is an unknown. Absent stagnation, this can be rectified by 
considering a change of variables (x, y) i— >• {q,p), where 

q:=x, p := -ip{x,y). 

This procedure is known variously as the semi-Lagrangian transformation or the 
Dubreil-Jacotin transformation. The effect is to map a single period of SI into a union 
of rectangles D = D^'uD' 2 ' C K 2 , since <9S! is mapped to the sets {p = 0}U{p = po}. 
Note that, by definition, {p — p{\ is the image of X under the transformation. The 
image of the air region Ox is thus 



Di := {(q,p) G D : < q < L, p 1 < p < 0}, 
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while the water region is mapped to 

D 2 ■= {(q,p) eD:0<q<L,p o <p< Pi}. 

With that in mind, we put 

T:={p = 0}, J:={p=pi}> B~{p = Po }. 

Let h — h(q,p) be the height above the bed of the point with x = q and lying on 
the streamline {ip = —p}, 

h(q,p) ■■= *)(q,p) +d, 

where rj = rj(q,p) is the vertical variable y in the new coordinates. More explicitly, it 
is the unique solution to 

~>P(q,t)(q,p)) = ~P, 

the existence of which is guaranteed by the absence of stagnation. By adopting the 
semi-Lagrangian coordinates, we are strongly exploiting the fact that there are no 
critical layers in the flow. 

Equation (2.2) can be reformulated as an equivalent problem for h: Find (h, Q) 
with h even in q, h p > in D, and satisfying the height equation 



(1 + h q )h pp + h qq hp - 2h q h p h p 



-h 3 p ~f(-p), in£>iUL> 2 , 



h = 0, 

h = £ + d(h), 
Here d is the depth operator 



on p = pi, 



onp=p , 
on p = 



(2.5) 



d(h) := 



L/2 



L/2 



Hq,Pi)dq- 



Note that we do not specify the value of d in advance, rather the correct value of d(h) 
emerges self-consistently from the equations. These equations can be found by applying 
the procedure as in [22] to obtain the interior equation, and using the following change 
of variables formulas to reformulate the jump condition: 



h q = 



d x = Og 



-d 



y/ghp' 
1 „ 



1 

y/ghp'' 



Note that the slope h q is continuous across the interface so could also be extracted 
from the jump condition on p = p\. The relative circulation can be calculated from h 
by 



r re i(p) 



L/2 x 

-L/2 



^dq, 



for pi < p < 0. 
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Lastly, let us set down some notation and describe the regularity of the solutions we 
wish to study. For any k G N, a € (0, 1), and smooth region 71 C M 2 , the space 
Cp+ a (7Z) is defined as the set of C k+a (7Z) that are L-periodic and even in their 
first coordinate. We are seeking smooth solutions to the problems enumerated above. 
Specifically, we look for solutions to the Euler problem of class y, the stream function 
problem of class y, and the height equation of class y" defined as follows 



(Q^, V )ey' := 

(Q,h) e y" :=: 



(c« r (n) n c£+°(fi\i)) 3 x c p 2 +«(M), 



i x (eg? (no n c^ a (n 2 ) n c p Q or (n) n c£?(n \ i)) 
x c 2 + a (R), 



per 

x(Cp 



jH , r (do n cli a (D 2 ) n c« r (D) n c 2 + a (D \ /)) 



Put more plainly, we want, e.g., solutions to the height equation to be C a in the 
whole domain, C 1+a up to the interface, and of class C 2+a away from the interface. 
The regularity of the other quantities is a direct consequence of that choice. 
Lemma 2.1 (Equivalence). The following statements are equivalent. 

(i) There exists a solution of class y to the steady stably stratified Eulerian 
problem (1.9)— (1.21) without stagnation. 

(ii) There exists a solution of class y' to the streamfunction problem (2.2). 

(iii) There exists a solution of class y" to the height equation problem (2.5). 
Proof. This lemma is routine. See, for example, [7, Lemma 2.1] or [22, Lemma 

2.1]. One point worth mentioning is that, while in our discussion of the formulation 
we only stated that tp was continuous across the boundary, as a solution of (2.2) it 
must in fact be of class Cp Cr (f2) by elliptic regularity (see Theorem A. 2). □ 

3. Local bifurcation for irrotational gravity waves. We begin by consider- 
ing the simplest case where the flow in both the air and water regions is irrotational 
(7 = 0). In this setting we need only consider the value of r re i on the lid. That is, 
the presence of wind here we interpret as having a nonzero r ro i on p = 0. Without 
loss, therefore, in this section we redefine r re i to be a positive constant. The height 
equation simplifies to: 



(1 



h 2 q )hpp 



h h 2 



*Zhqhphpq 



hi 



h = 0, 
h = £- 



+ 2g \p}h-Q = 0, 



d(h), 



in DiUD 2 , 



on p = pi , 



on P = Po, 
on p — 0, 



(3.1) 



and the relative circulation on the lid for a solution h is given by 

1 r L / 2 



L 



-L/2 



h p {q,Q)- x dq 



This follows from the fact that h q = on the lid. 



Our main theorem on this topic is the following. 

Theorem 3.1 (Local bifurcation for irrotational flows). Consider the existence 
of steady waves with a lidded atmosphere and irrotational flow in the air and water 
regions where F re / is given by 



\pi 



r, 



(3.2) 
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Assume that the following ideal local bifurcation condition holds 

-,gM + r I 2 cl coth(-^) <o. (ilbc) 

J- rel 

Then there exists a continuous curve of non-laminar solutions to the height equation 
for irrotational flow (3.1) 

C[^ = {(Q(s),h(s)):\s\<e}cy", 

for e > sufficiently small, such that (Q(0),/i(0)) = (Q(A*), H(X*)), and, in a suffi- 
ciently small neighborhood of (Q(X*), H(X*)) in J?", C[ oc comprises all non-laminar 
solutions. 

Remark 1. Using the equivalence of the three formulations, the theorem can be 
stated in terms of the original Eulerian problem as follows. Consider the existence 
of steady waves with a lidded atmosphere and irrotational flow in the air and water 
regions. Fix the (pseudo) volumetric mass flux in the air region to be p\, the height 
of the lid to be t; this forces any laminar flow to have (pseudo) relative circulation 
T re i defined by (3.2). //, in addition, the ideal local bifurcation condition (ILBC) is 
satisfied, then there is a corresponding continuous curve 

C\ oc = {(Q{s),u(s),v(s),g(s),P(s),i](s)) : \s\ < e} 

of small amplitude solutions to the Eulerian problem for an ideal fluid, which likewise 
captures all non-laminar solutions in a sufficiently small neighborhood of the point of 
bifurcation. 

It is also worth mentioning that, in contrast to the existence theory developed in 
hypotheses (3.2) and (ILBC) are both necessary and sufficient for local bifurcation 
to occur. 

The proof of this is result is developed over the next several subsections. 

3.1. Laminar flows. A laminar flow is one in which the free surface is un- 
perturbed, meaning that n = 0, and where the streamlines are parallel to the bed. 
In terms of the height equation formulation, this entails a solution with the ansatz 
h = H(p), and H(Q) = d(H). For such solutions, the PDE in (3.1) reduces to an 
ODE: 

H pp = 0, in (p ,Pi) U (pi,0), 
iH- 2 \+2g{p\H-Q = Q, on p = Pl , 
H(0)=i + d(H), 
H(p Q ) = 0. 

This problem can be easily solved explicitly, which leads to the following lemma. 

Lemma 3.2 (Laminar flow). For a fixed po,pi,£, and [p], if T vc \ is given by 
(3.2), then there exists a one-parameter family of solutions {(H(-; X),Q(X)) : X > 0} 
to the laminar flow equation (3.3) with H p > and where each solution has relative 
circulation F rc i on the lid. Explicitly, 



H(p;X) = { 



— r , J3l<_p<0 

1 rel A 

(3.4) 



P-Po 
X 



Pa < P < Pi 
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Q(A) = ^M(gi^M + r^-A 2 . (3.5) 

A 



In particular, the depth of the fluid at parameter value A is 



d(X) = — , (3.6) 
Pi - Po 

and the width of the corresponding channel is 

W(X) :=£ + d(X) = £+ — - — . (3.7) 
Pi - Po 

Remark 2. Let us make a few comments on this lemma. 

1. The proof of this lemma is straightforward. Notice that (3.2) is required to 
ensure that the resulting solution is continuous across the interface. Thus it is a 
necessary condition for the existence of laminar flows. 

2. Writing the corresponding solution in Eulerian form gives (u,v) = (U(y),0), 

where 

c-^L, 0<y<l 

U(y;X) = { 

z -d < y < 

\/ Pwater 

Thus the parameter X essentially dictates the relative speed in the water region while 
T re i dictates the relative speed in the air. 

3. Differentiating (3.5) in X, it is clear that X M> Q(A) is concave and has a 
unique maximum at X = Xq, where 

A§ = -5[p](Pi-Po). (3-8) 

3.2. Linearized problem. Fixing A > 0, we next linearize the full height equa- 
tion problem (3.1) around (H(-; A), Q(X)), which results in the following: 

m pp + m qq Hp = 0, inL>iUL>2, 

[iy- 3 m p ] =g\p\m, onp = pi, , g g . 

m = on p — po, 

m ~ d(m) =0 on p = 0. 

For simplicity, we specialize to L = 2ir; the general case can be approached via 
rescaling. Now, since we seek solutions that are 27r-periodic and even in q, it is natural 
to take m to have the ansatz m(q,p) — M{p) cos (nq), for some n > 0. Inserting this 
into (3.9) we see immediately that, for n > 1, M must satisfy the following: 

a\M pp = n 2 M in (po.Pi) U (pi, 0) 

{alMpj =g\p}M an p = Pl , , , 

M = on p = po, { n> 

M = on p = 0. 



Here we are denoting 

a\(p) ■= H p (p;X)~ l = 



Trcl Pi < P < 0, 
X p a < p < pi . 
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On the other hand, when n = 0, the equation in the interior are the same as above, 
but the depth operator d does not vanish, meaning that the boundary condition is 
nonlocal: 

a\M pp = in {p ,Pi) U (pi,0) 

blMp] = g ipj M on p = p u 

M = on p = po, { 0> 

M(0) = M(pO. 

Solving (P„) for any value of n > 1 will produce a 27r/n-periodic solution to 
the linearized problem. The next lemma gives the relation between the wavelength, 
circulation, and the parameter value A (which, recall, is associated with the speed in 
the water region). This can be seen as a form of dispersion relation. 

Lemma 3.3. For each n > 1, there exists a nontrivial solution M to (P n ) if md 
only if X = X n , where \ n satisfies 

= T rel coth (— ) - (AJ coth( ). (3.10) 

n 1 re i A* 

Such a A* will exist if and only if 

^i-r ?ol coth ( ^)<0. (3.11) 

n r rol 

Indeed, if it exists, X n is unique. If (3.10) (or, equivalently, (3.11) ^) holds, the space 
of solutions is one- dimensional and spanned by 

f sinh(^) Pl <p<0, 

^ nsmh( x , ) Po<P<Pi, 

where 

■ npt . 



sinh ( 



fi = n(X,n):= Trel . (3.13) 

sinh( " (Pl 7 P0) ) 
A 

Finally, there are nontrivial solutions to (Po) if and only if X = Aq ? where Xq is as in 
(3.8). ' 

Proof. Fix n > 1 and consider (P n )- By the ODE satisfied by M, we have 
immediately that 

M{p) = C\ exp (na^ p) + C*2 exp (— na^p), 

where 

Ci = ( C S in(Pl '° } ' for i= 1,2. 

[ C> m (po,pi), 

From the boundary conditions on the top and bottom, we have 



16 



S. WALSH, O. BUHLER, J. SHATAH 



Continuity of M at the interface implies that 

. , r npx ( 2 ) ,np . f n(pi-p ). 
C{ smh (- — ) = C{ exp (— — ) smh ( ). 

1 rcl A A 

Incorporating these observations, we can write M in the simplified form 

M«0) = Csinh(^-), 

; , (3.i4) 

M(^) = ^sinh( n( ^ l) ) 

where p, = /i(A,n) is as defined by (3.13). 

Lastly, we must ensure that the jump condition at the interface is met. We 
compute from the above expression for M that 



HM P j 



r r 3 el MW-A 3 Mf 
= CnT 2 rcX cosh(^-) 

1 rcl 

-/,CnA 2 co S h( n(Pl ~ Po) ). 

A 

Equating this to g [p] M(pi) found via (3.14) and simplifying yields 
9 M sinh (— i) = nr 2 cl cosh (— ^) 

-L rcl J- rc ] 

v2 , ,«(Pi-Po)« 



— /in A cosh ( ) . 

A 



Recalling the definition of u, this becomes 



M -r? cl coth(^)-A 2 coth(^f^l), 



n rcl v r re / v A 

which is the stated dispersion relation (3.10). 

Fix n > 1 and consider the map A e R + M> A 2 coth(n(pi —po)/X) G K + - Ele- 
mentary calculus confirms that it is a strictly increasing and nonnegative. Therefore, 
provided that (3.11) holds, i.e., 

^-r? cl coth(^)<o, 
n r re i 

there is a unique A = A* for which the dispersion relationship (3.10) is satisfied. On 
the other hand, if this inequality does not hold, then there will be no such A and thus 
no nontrivial solutions to the eigenvalue problem. 

Next consider the zero-mode case n = 0. Letting m(q,p) = M(p) in (3.9) we see 
that M must be piecewise linear. Moreover, the condition at p = implies 

Af(0) = d(M) = M{pi). 

From this we infer that M = M(0) on the interval [pi,0]. Due to the boundary 
condition at po and the piecewise linearity of M, we know 

M(p)=M(p 1 )^^, pe[ Po ,Pi}. 
Pi - Po 
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Using this to evaluate the jump condition reveals that 

A " -M(p 1 )=g\p}Mfa). 



Pi -Po 

Thus there is a nontrivial zero-mode solution if and only if 

This completes the proof. □ 

One final technical point needs to be made: since the laminar curve is parame- 
terized by A, and the solutions we seek depend on Q, we need to ensure that at the 
point of bifurcation Q is an invertible function of A. This is demonstrated in the next 
lemma. 

Lemma 3.4. For each n > 1 such that (3.11) holds, Q is an invertible function 
of A in a neighborhood of A* . 

Proof. Fix n > 1 satisfying (3.11) and let A* and M* be defined as in (3.10) and 
(3.12), respectively. From the formula (3.5), it is obvious that Q is a strictly concave 
function of A, and hence to prove the lemma it suffices to show that A* =^ Ao, where 
Ao is the unique critical point of Q. 

Multiplying (P n ) by M* and integrating by parts, we obtain 

° a 3 K {d p M* n f dp + n 2 f a K (M*) 2 dp + g [p] M* ( Pl ) 2 = 0. 
Upon regrouping terms, this becomes 

,2 / I HjT*\2 „ ir„H ii/T*^. \2 1 / „3 1 pi 7i / r*\2. 



0>-n 2 a K {M*y=glp\M*{ Vi y+ a 6 K (d p M* n ) 2 dp. (3.15) 

Jpo Jpo 

On the other hand, as M*(p ) = 0, 



\ /„ ( P 1 r ( / ^ (d p M* n ) dp) < ( r al o {dpM* n f dp) ( dp 
Po / \Jpo / \Jpo 

Pl al o {d p M:) 2 dp) 
po / V A o 



1 



Pl al o (dpM:) 2 dp. 



Po 

From this it follows that 

rpi 

< g M M^p,) 2 + / a\ {d p M* n ) 2 dp. (3.16) 

Inequalities (3.15) and (3.16) cannot be reconciled if Ao = A*, allowing us to conclude 
that this is never the case. □ 

3.3. Proof of local bifurcation. The objective of this section is to apply the 
theory of local bifurcation from simple eigenvalues to construct small amplitude (non- 
laminar) solutions to the wind wave problem, eventually culminating in Theorem 3.1. 
The machinery we employ is the classical work of Crandall and Rabinowitz, which, 
in the interest of readability, is included in the appendix as Theorem A.l. 
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Our first task is to put our problem into the framework of Theorem A.l. One 
cosmetic difference is that we wish to bifurcate from the family of laminar solutions, 
whereas Theorem A.l concerns bifurcation from solutions of the form (A, 0). With that 
in mind, let (h, Q) solve the height equation, and suppose h(q,p) = H(-; A) + m(q,p) 
and Q = Q{\). Then 



(1 + m 2 q )(m pp + H pp ) + m qq (m p + H p ) 2 
—2m q (H p + m p )m pq = 



in D 1 U D 2 , 



1 + m 



2 n 



m 



+ 2g {pi (m + H)-Q = 0, onp = Pl , 



m = 0, 

m + H -£- d(m) - d(H) = 0, 



on p = p Q , 
on p = 0. 



This can be restated cquivalcntly as 

F(\,m) = 0, 

where T = (J 7 !, J2, J3, J4) : K x X — > Y is defined by 

Tt(X, w) := (1 + (<>) 2 )(™« + H pp ) + + H p ) 2 

J- 2 (A, w) := (1 + (w&) 2 )(w$ + H pp ) + w$(wP + H p ) 2 
-2 Wq 2 \H p + Wp V)w$ 
l + w 2 



7" 3 (A, w) 
Fi(\w) 



2.9 M (w + H)\ I + Q 



(w — d(w) + H — d(H) -£)\ T - 
Here, the Banach spaces X and Y = Y\ x Y 2 x Y3 x Y4 are 

X := {h G C 2 p + a (D\I) n C« CI (D) : h( Po ) = 0, ft« G C^+ Q (A)}, 



(3.17) 



(3.18) 



Yi := C p 2 +«(^i \ I) n Cp er (Di), F 2 := C 2 p + a (D 2 \ I) n Cp 0r (D 2 ) 



^3 Cp CI .(7), F 4 :- C 2 + Q (T). 

Observe that, by Lemma 3.2, J-"(A,0) = for every positive A. In particular, we shall 
consider bifurcation from the lowest eigenvalue of the linearized problem found in the 
previous section. We shall therefore assume that (3.11) holds for n = 1 and denote 
A* := A*. 

For later reference, we now record the Frechet derivative of T with respect to w 
at (A*,0). 



F iw (\*,0)<p = (d 2 + H 2 2 ) tp® for i = l,2, 

= 2 [tf-Vp] - 2g M V, 
7" 4w (A*,0)(/? = (<p-d(tp)) T . 



(3.19) 
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Lemma 3.5 (Null space). The null space of J- W (X* ,0) is one-dimensional. 
Proof. Let <p G Af(J- w (X* ,0)) be given. By evenness, we can express (p via a 
cosine expansion: 

oo 

= ^2 (p n (p) cos (nq). 

n=0 

It follows that, 

J r [U (A*,0)(v3 n (p)cos(?ig)) = 0, n > 0. 

Equivalently, we must have that ip n solves (P n ). By Lemma 3.3 and the definition of 
A*, we know that tpi is nontrivial and that ip n vanishes identically for n^l. We have 
therefore shown that Af(T w (X*, 0)) is one-dimensional, and spanned by (p* := ip±, the 
unique solution to (P„) for n = 1. □ 

Now that we have ascertained the dimension of the null space, the natural next 
step in showing that T is Fredholm of index is to prove that the range is the 
(weighted) orthogonal complement of the null space. 

Remark 3. For the case we are studying where the air is irrotational, this is not 
particularly difficult if one takes the following approach. To study the solvability of 
J- W (X* ,0)ip = A for A € Y , we may project onto the individual modes by expanding 
PiQiP) = J2n Vnip) cos (nq). Doing so, we see that tp n solves an inhomogeneous 
version of the linearized problem (P n ). When n > 1, the problem can be viewed 
simply as two second- order ODEs with Dirichlet boundary conditions on the top and 
bottom, which must then be matched so that the jump condition on the interface is 
satisfied. The fact that the range has codimension 1 will arise as from this matching, 
and classical existence theory for ODEs. As before, the zero mode is must be dealt 
with using the fact for A* there is no eigenvalue. This procedure gives solutions 
of class C 2 away from the interface; the proof that the solution is in X follows from 
elliptic regularity, as described in Theorem A. 2. 

In §4> however, we allow the atmosphere to be rotational, and one of the effects 
of the vorticity is to make the linearized problem corresponding to (P n ) onerous to 
solve explicitly. With that in mind, it seems that a better adapted approach is to avoid 
separating variables, relying instad on purely PDE existence theory. Done this way, 
the proof of the range lemma is nearly identical in both regimes, and so this is the 
manner in which we have chosen to present it here. 

Even using PDE methods, though, the zero mode is somewhat special, since it is 
only there where the nonlocal operator d can be seen. For that reason, we will still 
wish to begin by projecting elements of X and Y onto their zero modes. We adopt 
the following notation, 

OW) == ^ r 9(q, ') dq, for any g E X, Y lt Y 2 , Y 3 ,Y 4 , or Y 5 . 

27T J -TV 

Lemma 3.6 (Range). A = (^1,^2,^3,^-4) € Y is in the range o/J" w (A*,0) if 
and only if it satisfies the following orthogonality condition 



if a 3 A!(p* dqdp+ [[ a 3 A 2 (p* dq dp + i / A S (f* dq+ [ a 3 A 4 (p* p dq = . (3.20) 

J J D-l J J D 2 2 J I JT 
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Proof. We begin by demonstrating necessity. Suppose that A £ 7£(.F W (A*, 0)). 
We may therefore let ip be given such that ^(A*, 0)<p = A. It follows that 



(cPtp* ,Ai) L2{Di) + (a?<p*,A 2 ) L2{D2) = II a 3 (tp P p + Hfaqq) ip* dq dp 

+ // a 3 (ifipp + Hpifiqg) ip* dqdp 

J J Do 



D. 

3 

u 

D 



a ippippdqdp 
// a 3 ip p ip p dqdp 

J J D 2 

a 3 (j)pip*\ dq 

a w* qq dq dp. 



(3.21) 



-Diun 2 



Here we have exploited periodicity and the fact that tp* vanishes identically on T. 
Continuing with the computation, 

(a 3 ^*,A 1 ) L2{ ) + (a 3 p*,A 2 ) L2{ } = // {a 3 p* + al p qq ) l pdqdp 

y ' ' JJd 1 ud 2 

h 3 <Pt<Pl - I^VpV*}) dq 



a A^ifp dq — d((p) / a ip*dq 

T JT 

glp}<p*<pdq 

f J (~Aa-gM <f\ V* dq 
a 3 AA<pp dq. 



Simplifying, we see that, indeed, the orthogonality relation (3.20) must be satisfied. 
The proof of necessity is complete. 

The next (more difficult step) is to show that (3.20) is sufficient. Suppose now that 
A satisfies (3.20), we wish to prove that there exists p £ X such that J 7 W (X*, 0)ip = A. 
First we consider the zero mode problem found by applying the projection P to the 
equation, which gives 



(3.22) 



Ai 


-dip 


for pi < p < 0, 


A 2 


= d 2 p <p 


for p a <p<pi, 


As 


= 2 [« 3 


Pp\-2glp\p, 


A 4 


= (¥- 


d{Tp)) T , 



where 

P~¥p, A = (A 1 ,A 2 ,A 3 ,Ai) :=0Mi, 
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Hence, 



¥ 1} (p) = 

?(2) 



rP r r 

/ / Ai{s)dsdr + dp + C 2 , 
Jpi Jpi 



(3.23) 



Tp^ip)= A2(s)dsdr + C 3 (p-p ), 

J Pa J Pa 

for some constants C\ , C 2 and C3 . The condition on T tells us that 

_ f° f p — 

- C lPl =A 4 - / Ax{r) drdp. (3.24) 
Jpi Jpi 

Continuity across the interface implies 

C lPl + C 2 - ( Pl - Po )C 3 = [ 1 [ A(r) dr dp. (3.25) 

Lastly, the jump condition equates requires that 

rf el d - (A*) 3 ^ A 2 (p) dp + C^j-g \ P \ (C lP + C7 2 ) = pt 3 . (3.26) 

Collecting together (3.24), (3.25), and (3.26), an easy calculation reveals that unique 
solvability of the zero mode problem is equivalent to 

-3M(ft-ft)^(A*) 3 . 

But this is just the statement A* 7^ Ao, which is proved in Lemma 3.4. We have 
therefore shown that the zero-mode problem has a unique solution. 

Returning to the question of solvability of the full problem, we observe that, in 
light of the previous analysis, it suffices to assume 

A e Y Q := (1 - F)Y = {B e Y : FB = 0}, 

and to solve in the space 

ip e X := (1 - F)X = {(/)eX :F(j) = 0}. 

In fact, this means that we may take d((p) = 0. 

We shall approach the question of solvability incrementally. Fix e > 0, and define 
: X Q -> Y by 

= (e- T lw (\*,0), e- T 2w {\*,0), - JF 3 «, (A* , 0) , -T iw (\* , 0)). 

First consider the approximate problem: 

C^^^A. (3.27) 

Claim 1. For a sequence of e > tending to 0, there exists a unique solution ip^ 
to (3.27). To see this, note that by a method of continuity argument, the solvability 
of (3.27) is equivalent to that of the equation 

£0V £) = A, 
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where 



C^i :=C[ € \ forz^3, £W 3 </? (e) := 2 



That is, we may safely ignore the zero-th order term on the interface. Let £ 6 
Cp er (D)nCp+ Q (Di)nCp+ a (Z?2) be any function that exhibits the following properties: 

e| fl = 0, {a%] - -A3, £|t - -A*, and P£ = 0. 
Then, the solvability of (3.27) is equivalent to that of 

£W^ (£) = A, (3.28) 

where ip^ :— ip^ — £, and 

A:= (A -/ft M-C^k, 0, 0). 

The Fredholm solvability of (3.28) follows from linear elliptic theory (cf. Theorem 
A. 2); we therefore need only to prove that the homogeneous problem associated with 

£( £ ) has no nontrivial solutions. Suppose that <j> solves the homogeneous problem for 
some fixed e > 0. By evenness and the fact that ¥(f> = 0, we can expand <f> in a cosine 
series of the form 

oo 

<t>{<l,p) = ^2 <t>n(p) cos (rig). 

71=1 

Applying the operator £( e ) to <j> reveals that tj> n satisfies the eigenvalue problem 
dl<j) n = (e + n 2 H^ n , [a 3 <Vn] = Mpo) = <f> n (0) = 0. 

This is a consequence of the fact that d(cf> n ) = 0, as n > 1. This equation can be 
easily solved explicitly, and one can readily see that, for generic e small, there are no 
nontrivial solutions for any n > 1. We omit the details in the interest of brevity. The 
first claim is proved. 

We may therefore let e„ be a positive sequence converging to as in Claim 1, and 
consider the corresponding sequence of solutions to (3.27), call them {tp^}. 

Claim 2. {ip (n) } is bounded uniformly in C*+ a (D?) n Cp-+ Q (I^). We argue by 
contradiction. Suppose that {<p( n >} is not bounded uniformly. Possibly by passing to a 
subsequence, we may suppose that ||^' n ^|lc l +«(i^ + ||V'^llc^ (i5^) — ^ 00 • bet := 
^/(II^Hll^^ + II^Hll^^). Thus has unTt C^ a (D^nC^+ a (D^- 

norm, and, by linearity, is a solution to 

£( e »U(«) = =■ A n) 

IM n) llc^(Dl) + 11^11^(153) 

Since A") converge to in Cp+ a (-Di) fl C^ r a {D 2 ), we may extract a subsequence 
{</)("*)} converging to e C 2 (D\I) n C 1 ^) n C 1 ^), a classical solution of 

^(A*, 0)0 = 0. 
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This is achieved by using Schauder-type estimates for the approximate problem, and 
the compactness of the embedding of C k+a into C k on bounded domains. From 
Lemma 3.5 it follows that = vip*, for some i/gR, Finally, we note that it must be 
the case that H^Hci+ajT^j + ll ( / ) llc 1+cv (i>j") ~ ^> an< ^ nence v 0- 
Recall that, by definition the of we have 

aV£ (e "V n) = a 3 <p*A. 

Integrating over D, we obtain 

(aV,-4i) L2(Dl) + («V,^ 2 ) L2(D2) = // aVMW-^-ff,t^* 



D 2 



a 3 <p*(e n <p^-^-H 2 p ^)dqdp 



aVV (n) dqdp 



D1UD2 



r/<7 



DiU-D 2 



— - y a ^393* dq — J a Anp* p dq. 

By the orthogonality condition (3.20), all the terms involving A cancel, leaving only 
the statement that 



a"W nJ dqdp = for all n > 1. 

'D1UD2 

As an immediate consequence, we have 

a 3 <f*(j) dq dp = 0, 

'DiUfl 2 

which contradicts the fact that 4> — VL P* ■ This completes the proof of the second claim. 

Up to this point we have shown that there exists solutions {^"'j to (3.27) for a 
sequence of e„ — > that are uniformly bounded in C l+a . By passing to a subsequence, 
we have that there exists a weak solution to the original problem J r w (X* ,0)ip = A. 
Elliptic regularity ensures that the weak solution is actually strong, and, in particular, 
is an element of X (cf. Theorem A. 2). It follows that A is in the range. The proof of 
the lemma is complete. □ 

Lemma 3.7 (Transversality). The following transversality condition holds: 



Here ip* denotes a generator of A/"(J r u ,(A*, 0)). 



(3.29) 
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Proof. In view of Lemma 3.6, it suffices to show that A := Fxwi^*, 0)<p* fails to 
satisfy the orthogonality condition (3.20). That is, if we put 

,3 /!_,„* j„j„ 1 //(.,„* j„ / „3 



a e?g<ip + yy aA2(p*dqdp+- j A^tf* dq + j a A^Ppdq, 

we must prove S ^ 0. To do this, we first compute 
•7 7 iAt«(A*,0)</3* = 0, 



2 

(A* 

7- 3 a w (A*,0)^ = (-3(A*) 2 (^)^) jI 



^4At«(A*,0)^* = 0. 

By the equation satisfied by (p* , we see that 
Thus, 



(A*) 2 (^V;) p = (A*) 2 (^) 2 + (^) 2 , in£> 2 . 



From this we deduce, 



a a A 2 (f*dqdp = 2 // (ip*y dqdp, (3.30) 



and 



A, f : ^ = -M^l! p* dq 



1 

2.1 ' 2 



(^p) dqdp- - // (</) dqdp. 
D 2 ^ •/ 



(3.31) 



In light of (3.30)-(3.31), and the calculated values of the remaining components of .4, 
3 = ^/X ^fdqdp-l(Y) 2 JJ d {V* p ?dqdp. (3.32) 
Recall that {ip*)^ takes the form 

y*(p) =Msinh( ° ) in D 2 , 
for an explicit constant /i. Hence, 



and thus 



(A*)V» 2 = M 2 cosh 2 (*^° ) in D 2) 



<p*(p) 2 - 3(A*)V» 2 = - M 2 (l + cosh 2 (^)) in ^2 
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From this we conclude S < 0. □ 

With these lemmas, Theorem 3.1 becomes a simple consequence of the Crandall- 
Rabinowitz bifurcation theorem. 

Proof of Theorem 3.1. Assuming (3.2) and (ILBC), we are justified in defining 
J- and A* as in (3.18) by Lemma3.2 and Lemma 3.3, respectively. To complete the 
proof, we need only confirm that the hypotheses of Theorem A.l are satisfied. But 
parts (i) and (ii) of the Crandall-Rabinowitz theorem are obviously true. Lemma 3.5 
and Lemma 3.6 together give (iii), while (iv) was proved in Lemma 3.7. □ 

4. Local bifurcation with vorticity in the atmosphere. Next consider the 
case where the density remains constant in each region, but we assume that the 
air region has a nontrivial vorticity strength function. The height equation for this 
scenario is 

(1 + h 2 q )h pp + h qq h 2 p - 2h q h p h pq = -^{-p)hl inD 1 \JD 2l 

+ 2g [pj ft. - Q = 0, onp = pi, , 41 , 

h = 0, onp = p , 

h = £ + d(h), onp = 0. 

Written in the new variables, the relative circulation becomes 

1 [ L / 2 1 + ft, 2 
r rc i = y — r — - dq, for p 1 < p < 0. 

■ Ij J —L/2 Up 

Our main theorem is the following. 

Theorem 4.1 (Local bifurcation with atmospheric vorticity). For given pi, I, 
and 7, define F rc i by 

d p {TUpf ) = 27(-p), l = J° f^p-y ( 4 - 2 ) 

If the local bifurcation condition holds (cf. Definition 4-5) then there exists a contin- 
uous curve of non-laminar solutions to the height equation for irrotational flow in the 
water and vorticity strength function 7 in the air (4.1) 

C? oc = {(Q(s),h(s)) ER x X :\s\<e}, 

fore > sufficiently small, such that (Q(0), h(0)) = (Q(X*), H(X*)), a laminar solu- 
tion with (pseudo) relative circulation T re \, and, in a sufficiently small neighborhood 
°f (Q(A*), H (A*)) inRx X, C' loc comprises all non-laminar solutions. 

Remark 4. As in Theorem 3.1, we can interpret the above statement in terms of 
the Eulerian formulation. The resulting statement is as follows. Fix the period to be 
27r, (pseudo) volumetric mass flux in the air region pi, lid height t, and with vorticity 
strength function 7 and define T rc \ by (4.2). If the local bifurcation condition holds, 
then there is a corresponding continuous curve 

Cioc = {(Q(s),u(s),v{s), Q (s),P(s),r)(s)) : |s| < e} 



1 + 



26 



S. WALSH, O. BUHLER, J. SHATAH 



of small amplitude solutions to the Eulerian problem with an irrotational water region 
and a vorticity strength function 7 in the air region, which likewise captures all non- 
laminar solutions in a sufficiently small neighborhood of the point of bifurcation. 

Lemma 4-6 provides an explicit size condition (4.9) under which the local bifurca- 
tion condition holds. We have left it in the more abstract form here in order to give 
the most general statement. On the other hand, the compatibility condition is, in fact, 
necessary for the existence of laminar solutions. 

4.1. Laminar flows. Consider the laminar flow problem where we seek solution 
to (4.1) with the ansatz h(q,p) — H(p). Then the PDE simplifies to the following 



H pp — 

H.„ 



L pp 

2g\p}H-Q 
H(0) 
H(pa) 



in (p ,Pi) 
~f(- P )Hl in(pi,0) 



0, 



on p - 
d(H), 



Pi, 



(4.3) 



Evaluating the relative circulation gives, 

H p 1 = r rcl , 



in \pi,0] 



Again, this is explicitly solvable, but there are some compatibility conditions that 
are necessary to ensure continuity across the interface. 

Lemma 4.2 (Laminar lemma). If the compatibility condition (4.2) is satisfied, 
then there exists a one-parameter family of solutions {(if(-; A), Q(X)) : A > 0} to the 
laminar flow equation (4.3) with H p > and where each member of the family has 
prescribed relative circulation T rc \ . It has the explicit form 



H(p; A) 



clr 



pi 



r re i(r) 



Pi -Po 
A : 



P-Po 
A : 



Pi < p < 



Po < P < Pi 



(4.4) 



and 



Q(A) 



2.9 [ 



(Pi - Po) 
A 



r re i(pi 



A 2 . 



The depth of the fluid at parameter value A is 

d(X) = 

Pi - Po 

and the width of the corresponding channel is 



(4.5) 



(4.6) 



W{\) :=£ + d(\)=£ + 



Pi - Po 



(4.7) 



Proof. Examining the ODE in the air region reveals that 

f 



0, 



1, 



for p € (pi,0). 
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This implies directly that (T 2 cl )' = 27. Integrating the equation H 1 = r rc i, we have 



for some constants C. Since H pp vanishes in the water region, and H(p ) = 0, 



P-Po 

x • 



for some constant A > 0. Continuity at the interface then implies that C = (pi— P o)/X. 
Using this, and the fact that d(H) = H(pi), the condition on T becomes 



1 = 



dr 



Pl r re i(r)' 



which is the second part of the compatibility condition (4.2). 

Lastly, we use the jump condition on the interface to determine Q as a function 
of A: 

Q = [tf" 2 ] + 2.9 Ipj Hip,) = r rel (p x ) 2 - A 2 + 2. 9 Ipj , 

which is (4.5). □ 

Remark 5. The dependence of Q on X is essentially unchanged from the ideal 
case (3.5). In particular, it is concave with a unique maximum occurring at Xq as 
defined in (3.8). 

4.2. Linearized problem. Proceeding as before, we seek to linearize the full 
height equation around the curve of laminar flows. For a fixed A > 0, this gives 



2- 



(a 3 m p ) p + (am q ) q = 0, in D\ U D 

[a 3 m p | = .9 Ipj m, on p = pi, 

m — 0, on p = po, 

m — d{m) =0, on p = 0. 



(4.8) 



Here we are again using the convention 

a = a{p:X)=H p {p-X)-' = { \^ Pl < P < °' 
v ; pV ' \ A, po <p <pi. 

Note that the linearization of the circulation identity is 

1 r 

— / am p dq = for pi < p < 0, 

and hence each of the solutions of the linearized problem will have no relative circu- 
lation in the air region. 

We seek solutions with the ansatz m(q,p) = M(p) cos (nq). First consider the 
case where n — 0, i.e., there is no 9-dependence. Then 

(a 3 M p ) p = 0, in D t 

M pp = 0, in D 2 

[a 3 M p J = g \p\ M, onp = Pl , 

M = 0, on p = po, 

M - d{M) = 0, on p = 0. 
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Since M is linear in D2, the boundary condition at the bottom implies that it takes 
the form 

M^( P ) = P^Jto-Mfa), pe\p ,Pi}. 
Pi - Po 

On the other hand, since M(0) = d(M) — M(pi), we must have that M p vanishes at 
least once in the interior of D\. Since o?M v is constant in D\, we conclude that M^ 1 
must, in fact, be identically zero. Thus the jump boundary condition states 

— - — %i) = jM%). 

P\ -Po 

Just as in the irrotational case, we see that there can be a zero-mode solution if and 
only if A 3 = —g [p] (pi — po), that is, A = A (cf. Lemma 3.3). 
If we take n > 1, then (4.8) becomes 

(a 3 M p ) p = n 2 aM in (po,Z»i) U (pi, 0) 

K tf J = 5 M M onp = Pl , 

M = onp = j3 , " 

M = on p = 0. 

Notice that periodicity implies that e?(M) = 0. 

We shall approach the problem of finding solutions to (P„) using a variational 
method. Define the Rayleigh quotient & by 

A) := -0 — , A > 0, ip e < 

Jpo d P 

where the admissible set 

^ := {0 e i 2 (N,0]) n ^(bo,^)) n ^((pi.O]) : 0(0) = 0( Pl ) = 0}. 

A straightforward argument then gives the following lemma. 

Lemma 4.3. If, for fixed A > 0, ip is a critical point of&(-; A) and Si{jp\ A) = —n 2 , 
for some n > 1, then ip solves (P„) for this value of n. 

Using very basic estimates, we can show that for A sufficiently large, min^ A) 
will be greater than — 1. More generally, we can show that the minimum goes to —00 
as A — > +00. This is the content of the next lemma. 

Lemma 4.4. Let a m i n he the minimum value of a on [p\, 0] (which does not depend 
on X). Then, for each n > 1, if 

A 2 > a 2 mm - ~9 M , 

we have &(ip; A) > — n 2 , for every tp e sz# . 

Proof. Fix any A as in the hypothesis and let tp G srf be given. Then 

p0 
(a 3 Lp 2 + n 2 ap 2 ) dp > a min / ((a min ip p ) 2 + (nip) 2 ) dp 

1 Jpi 

> 2na min / pp<pdp = -2na 2 nin ^(pi) 2 . 
Jpi 
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On the other hand, since oA 2 ' = A, 



(a 3 (fil + n 2 aip 2 ) dp > 2n\ 2 ip(p 1 ) 2 



pi 

(,.•1 : ^\ 

'Po 

Summing these together and recalling how we chose A, we find 

(aVp + »W) dp > (2n\ 2 - 2na 2 min ) V >{p 1 ) 2 > -g [p] p(pi) 2 ■ 

Rearranging terms, this implies that M(<p; A) > — n 2 . □ 
Let us define 

^(A) := min M(ip;X). 

Then, if we can show that ^(A*) = — n 2 , for some A*, Lemma 4.3 guarantees the 
existence of a solution to {P' n ), and thereby a 27r/n-periodic solution to the linearized 
problem. We are most interested in showing that —1 is in the range of v. The 
preceding lemma provides a lower bound for v when A is sufficiently large. In order to 
guarantee that —1 is in the range of v, therefore, we need only verify that z^(A) < —1 
for some positive A. This will not be true in general, and so we are forced to make it 
a hypothesis. 

Definition 4.5. We say that the local bifurcation condition is satisfied provided 

that 

inf v(\) < -1, (LBC) 

A>0 

or, equivalently, if 

there exists a nontrivial solution to the linearized problem {P' n ) for n = 1. (LBC) 



These are abstract conditions that are both necessary and sufficient for local 
bifurcation. One way to derive an explicit sufficient condition is to require that 
£i$((p; e) < — 1, for some e > and a conveniently chosen cp E . This is the approach 
of the next lemma. 

Lemma 4.6 (Size condition). Suppose that the prescribed circulation T re i, the 
pseudo-volumetric mass flux in the air region pi, and the density jump [p] collectively 
satisfy the following size condition: 

9 M Pi + / (r ro i(p) 3 + p 2 T Ici (p)) dp < 0. (4.9) 
Jpi 

Then (LBC) holds. 

Proof. This can be seen easily by taking 



<P(P) ■-- 



P/Pi, Pi <P < 

(p-Pi)/(Pi -Po), Po<P<Pi- 



and evaluating lim^o &(<p; A). Here we have made use of the fact that a^ 1 ' — L ro i to 
express the size condition only in terms of prescribed quantities. □ 
Lemma 4.7 (Monotonicity of v). If v(\) < 0, then dv{\)/d\ > 0. 
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Proof. Let Cip := — (a 3 ip p ) p . Then for each A, let tp £ &/ solve the problem 

Cip = v(X)aip, in (p ,Pi) U (p!,0) 
[a 3 ^p] = g Ip} <p, on 7, 
= onTUB. 

From this we compute 

dp = (3a 2 aip p ) p + vaip + vaip + vaip, in (po,pi) U (pi,0) 
l a3( Ppl + l3a 2 atpj = g [p] 0, on 7, 

9? = 0, on TUB, 

where a dot denotes differentiation with respect to A. Letting (•,•) be the standard 
L 2 -inner product, we have 

(dp, ip) - (dp, (p) = j>(a</3, tp) + v{aip, tp) + ((3a 2 dipp) p , <p) . 

Upon integrating the last term by parts, we find ultimately that 

(Cip, <p) - (dp, ip) = u(acp, tp) + v(aip, ip) - Z(a 2 a<pp, ip p ) - [3a 2 dv3 p ] tp(pi). 

However, if we simply compute the difference using integration by parts, we see 

that 

(C<p,<p) ~ (Ccp,(p) = {a 3 ip p j <p(pi) - [a 3 (^ p ] ipipi). 

Equating the last two lines and exploiting the jump conditions satisfied by ip and ip, 
we obtain the following identity 

v{aip, ip) + v(aip, ip) — 3(a aip p , (p p ). 

Observe that a > and a > (in fact, cS 2 ^ = 1 and — 0), hence the right-hand 
side above is positive. If v is negative, it must then be the case that v > 0. □ 

Lemma 4.8 (Existence of the minimizer). Suppose that the size condition (4.9) 
holds. Then there exists a unique value A* > such that v(\*) = — 1. Equivalently, 
there is a unique value of A for which there is a nontrivial solution to (4.8) with 
the ansatz m(q,p) — M (p) cos q. Moreover, Q is an invertible function of X in a 
neighborhood of A* . 

Proof. From Lemma 4.4, we have that v(X) > — 1 for A sufficiently large, while 
from Lemma 4.6 we know that v(\) < —1 for A sufficiently small. By continuity, 
there exists A* such that ^(A*) = — 1. Moreover, the monotonicity of v established in 
Lemma 4.7 implies that A* is unique. 

Next, since Q is a concave function of A according to (4.5), it will be a bijection 
locally so long as A* does not coincide with the unique critical point Ao of Q. Recall 
that Ao satisfies 

= -9M (Pi ~Po), 

or, put more suggestively, 

J I a A - 3 * = A -(p 1 -p o) = -J s . 
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Then, for any ip <G si, we have the estimate 

2 



/ rpi \ * / rpi \ / rpi 

^(Pi? = / 'Ppdp) < / a\y p dp\ / a^dp 

\Jpo / \>^Po / Wpo 



Po 

1 



9M 



a\y p dp- 



It follows that 



rpi rO 

< 9 M <p(pi) 2 + / a\y p dp<g Ipj <p( Pl ) 2 + / a 3 x y p dp. 

J pa Jpo 

As this implies 8%{y>; Ao) > for all ip <G si, we conclude that /z(Ao) 7^ —1, and hence 
A* 7^ A . This completes the proof. □ 

4.3. Proof of local bifurcation. Let us set h(q,p) = H(p;X) + m(q,p) and 
Q = Q(A). The equation satisfied by m is thus 

J*(A,m) = 

where T = (Ji, J" 2 , J3, J4) : K x A" ->• Y is defined by 

Ji(A, w) := (1 + «)) 2 )( W « + H pp ) + fflW^) + # p ) 2 

- 2™«(iJ p + + 7 (ff p + «;W) 3 

J- 2 (A, w) := (1 + («;< 2) ) 2 )(«$ + H pp ) + W ( 2 )( Wp 2 ) + H p f 

2w^(H p + w£))w$ + 7 (tf p + w^f (4.10) 



J* 3 (A, w) := - 



-2sM (w + H)U + Q 



(H p + w p ) 2 
J" 4 (A, w) := (to - + - d(JT) - | T . 

Here, the Banach spaces X and Y = Y\ x Y2 x Y3 x Y4 are 

X := {/* e C 2 +«(I5\ /) n C«(Z>) : h(p ) = 0, /»« e Ci+ Q (A)}, 



Yi := C 2 + a (£»! \ /) n c^pi), r 2 := c 2 + r a (£> 2 \ /) n c« r (D 2 ) 

Y 3 '■= Cper(-O) *4 : = Cptr a ( T )- 

For later reference, we now record the Frechet derivative of T with respect to w 
at (A*,0). 

T iw (A* , 0)p = (9 2 + # 2 9 2 + 3 7 ff 2 d p ) for i = 1,2, 

J- 3w (A*,0)^ = 2 [tf-Vp] - 2.9 M ip, (4.11) 
^4«>(A*,0)v3 = - rf(v)) T ■ 



Lemma 4.9 (Nullspace). T/ie null space of T w (A* , 0) is one- dimensional. 
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Proof. Let <p be an element of the nullspace. By the evenness built into the 
definition of X, we may expand ip via a cosine series: 



<P(Q,P) = ^2 fn(p) cos (nq). 



71=0 

ft follows that, 



7" w (A*,0)(y>„(p)cos(ng)) = 0, n > 0. 

Equivalently, we must have that ip n solves (-P,' ). By Lemma 4.8 and the definition of 
A*, we know that ipi is nontrivial. We have already seen that there are no nontrivial 
solutions of (-P,',) with n = 0, and hence ipo must vanish identically. If ip n ^ for some 
n > 1, then it belongs to the admissible set srf and hence &(p n ;\*) = —n 2 < — 1. 
This contradicts the definition of A* as the minimizer. Thus all the ip n with n > 1 
vanish identically. We conclude that the nullspace is generated by ipi = ip* . □ 

Lemma 4.10 (Range). A = (Ai,A2,As,A4,Ab) EY is in the range ofF w (\*,0) 
if and only if it satisfies the following orthogonality condition 

I a 3 A!ip* dqdp+ [ a 3 A 2 ip* dq dp + - / A 3 (p*dq+ / a 3 A 4 ip*dq^0. (4.12) 

J J D-l J J D 2 ^ J I JT 

Proof. First assume that A is in the range of T w (\* , 0). Then there exists ip E X 
such that J" w (A*,0)y3 = A. Writing ./^(A*, 0) for i = 1,2 in self-adjoint form, we 
have 



3 8 p (a 3 d p p {l) ^ + a- 2 d 2 q p (t) = A {l K i= 1 . 2. 

Thus 



(a 3 <p*, Ai) L2 ( Di) + (a 3 ip*, A 2 ) L2{L>2) = // ((a 3 tp p ) p + a<p qq ) if* dqdp 



((a 3 ip p ) p + a<p qq ) p* dqdp 
J J a 3 p p p* p dqdp 

-If a 3 p p p* p dqdp- / [a 3 ^*! dq 
JJd 2 ' J i 

atpip* qq dq dp. 



DiUD 2 

Here we have exploited the periodicity and the fact that ip* vanishes identically on 
T. This is precisely what we have in (3.21), and, since the conditions on the interface 
are unchanged, proceeding as in that proof shows that (4.12) is necessary. 

The sufficiency of (4.12) follows from a straightforward generalization of the ar- 
gument given in Lemma 3.6 which we omit. □ 

Lemma 4.11 (Transversality) . The following transversality condition holds: 

T Xw (X*,0)p* <?K{T w {\*,0)). (4.13) 
Here p* denotes a generator of M{T W {\* 



STEADY WATER WAVES WITH WIND 



33 



Proof. In view of Lemma 3.6, it suffices to show that A := Fxwi^*, 0)<p* fails to 
satisfy the orthogonality condition (4.12). That is, if we put 

S := I a 3 Aip* dqdp+ [ a 3 A 2 (p* dq dp + - / A 3 p* dq + / a S A±y*dq, 

J J Dt J J D 2 ^ J I JT 

we must prove 5^0. To do this, we first compute 
Ti\ w {\* ,<d)p* = 0, 



■ ?7 2A« 1 ,(A*,0)^* = - 7 — -(p^-— -<p*, 



2 * 67 
(Vp^w ~ (A* 

^(a*,ov* = (-3(a*) 2 (^;)( 2 )) j 



T iXw {X*,0)p* = 0. 
By the equation satisfied by y>* (P r ' ), we see that 

-dgip* = tp* = a~ 1 d p (a 3 dptp*). 

Thus, 

(A*) 2 (^V;), = (A*) 2 (^) 2 + (^) 2 , \nD 2 . 
From this we deduce, 

/ / a 3 A2P* dq dp — 2 / / {tp* ) 2 dq dp — 67 / / (p*<p* p dqdp 
J J D 2 J J d 2 J J D 2 

= 2 ff {p*) 2 dqdp, 



(4.14) 



and 



\fA^dq = - 2 



3 -^Jy P ) {2 Vdq 



(^;) 2 rf 9 dp-- // {p*ydqd P . 

z .1.ID2 ^JJD 2 



(4.15) 



In light of (4.14)-(4.15), and the calculated values of the remaining components of .A, 



1 / / , n 2 j j ^/\*n2 / / / *\2 



, (^y*Z*-~(A*r // HYdqdp. 

L JJD 2 z JJfl 2 

Since the equation satisfied by (y>*)( 2 ) is the same as for the irrotational atmosphere 
case treated in the previous section, we know that S < by the same argument as 
Lemma 3.7 (cf. (3.32)). This completes the proof of the lemma. □ 

Proof of Theorem J^.l. Under the assumption that the compatibility condition 
(4.2) and the local bifurcation condition (LBC) hold, Lemmas 4.9-4.11 verify that 
the hypotheses for Theorem A.l are satisfied. The conclusions of Theorem 4.1 follow 
immediately. □ 
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5. Local bifurcation with an unbounded irrotational atmosphere. In 

this section and the next we consider the situation where the atmosphere region 
has infinite vertical extent. We begin by supposing that the flow in both the water 
and air are irrotational. In the moving frame, the governing equations for the stream 
function formulation are as in (2.2) (taking pW to be constants, d = 1, and 7 = 0), 
with the only difference manifesting in the boundary condition at infinity: 

Aip = 0, in n, 

[I Wf] + 2 5 Ipj ( V + 1)-Q = 0, on y = V (x), 

ij) = 0, on y = r](x), (5-1) 

^ = ~Po, on y = -1, 

V ip — > (—A, 0) uniformly in x as y — > 00. 

Here A = \J p( 1 )(C7 — c) is the speed of the undisturbed wind in the co-moving frame 
multiplied by the square root o the density of the air. 

The change in the domain necessitates a change in the spaces in which we seek 
solutions. In particular, we must specify the behavior as y — > 00. As can be seen in 

(5.1) , for the irrotational regime we are interested the case where (u,v) — > (U,0), for 
some constant U < c as y — >• 00, which is equivalent to requiring that V^ip limits to 
(—A, 0), for some positive constant A. With that in mind, we define 

^0 := {(it, v, g, if) £ : 3U < c, (u, v) — > (U, 0) uniformly as y —> 00} 
^0 : = {(Q,^,V) e y' : 3A > 0, V- 1 ^ -> (-A,0) uniformly as y -> 00}. 

Theorem 5.1 (Local bifurcation with unbounded irrotational atmosphere) . Sup- 
pose that the volumetric mass flux in the water region po and the density jump [p] 
satisfy the local bifurcation condition 

pgcothl + fl[p] >0, (5.2) 

then the following statements are true. 

(i) There exists a continuous curve of non-laminar solutions to the stream func- 
tion equation for irrotational flow in the air and water, and unbounded atmosphere 
region (5.1) 

C[ oc = {(Q(s)^(s),r,(s))&^--\s\ <e}, 

for e > sufficiently small, such that (Q(0), i/>(0)) = (Q*, Wq), and, in a sufficiently 
small neighborhood of (Q*,^) in K x X, C[ oc comprises all non-laminar solutions. 
In particular, we have 

r](s)(-) = s cos (O+oOO inCl+ a (^)- 

(ii) There is a corresponding continuous curve 

C loc = {(Q(s),u(s),v(s),g{s),P{s),r,(s)) G ,y ■ N < e} 

of small amplitude solutions to the Eulerian problem which likewise captures all non- 
laminar solutions in a sufficiently small neighborhood of the point of bifurcation. 
Unlike to the explicit size condition (4.9) obtained in the rotational lidded regime, 

(5.2) is both necessary and sufficient. We caution, however, that (5.2) is stated in 
dimensional variables, and so some additional work must be done before attempting 
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to draw physical conclusions. For instance, note that we are taking d = 1, effectively 
setting the length scale of the system to match the depth of the ocean (which is 
typically greater than 3000 meters), and constructing minimally 27r-periodic waves. 
Clearly such waves are not typical. 

A simple generalization of Theorem 5.1 can be made where we allow an arbitrary 
depth d > 0, and perturbations of the form 2ir/k, k > 1. Let k* be the smallest 
natural number such that 

plk* coth (k*d) + gd 2 [p] > 0. (5.3) 

Then the conclusions of Theorem 5.1 hold, with the only modification being that the 
non-laminar solutions are of the form 

rj(s)(x) = s cos (k*x) + o(s). 

and the depth of the water region is d. Notice that k* — 1 and d = 1 if and only if 
(5.2) holds. 

5.1. The transformed problem. At a mathematical level, the major differ- 
ence in going from the lidded and to the unbounded case is the loss of compactness 
in the domain. If we were to attempt to use semi-Lagrangian coordinates and the 
height equation formulation, we would not expect that the operator T to be Fred- 
holm. This is a substantial technical obstacle, but one that has been grappled with 
extensively in the literature of traveling waves in oceans of infinite depth (cf., e.g., 
[12, 13]) and solitary waves (cf., e.g., [4]). Typically, the strategy is to rely on con- 
centration compactness arguments or a Nash-Moser iteration scheme. But aside from 
the inherent complexities of these tools, were we to attempt this approach, we would 
still be forced to require the absence of stagnation points and critical layers in order to 
justify the Dubreil-Jacotin transformation. This an especially restrictive assumption 
in the unbounded atmosphere regime. To see why, consider flows with a constant 
nonzero vorticity in the air region, which is the subject of the next section. A laminar 
flow ([/, 0) of this type will necessarily satisfy \U\ — > oo as y — > oo. Thus, depending 
on the sign of the vorticity, critical layers are expected, though they are neutered (see 
the remark following Lemma 6.2). 

Rather than adopting semi-Lagrangian coordinates, therefore, we shall employ 
a more robust (though less elegant) transformation to fix the domain. Let f^o = 

n { 1] U where 



fiW = [-7r,7r] x [0,oo), n {2) = [-w,ir] x [-1,0]. 

Suppose i](e, •) is a one-parameter family of free surfaces with 77(0, •) = 0. Then we 
may take r2(e) and tp(e, •) to be the corresponding families of fluid domains and stream 
functions, respectively. Let T = T(e,x,y) : il(e) — > f2o be a smooth diffcomorphism 
mapping Sl(e) to fio for each e > such that 

T(0, •) = tn (the identity map on £!o) 
T(e,x,j](x)) = (x,0) 
T(e,x,-1) = -1 
[T(e, •) — iji(e)] — >■ as y — > 00, uniformly in x and e. 
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Any such map will serve the purpose of fixing the domain, but for simplicity we 
make a particular choice. Suppose that in a neighborhood of the x-axis, T is just the 
flattening map: 

Ti(e,x,y) := ni(x,y) := x, 

T 2 (e, x, y) := V ^ L ^ x(e, y) + y(l - x(e, y)), 
1 + 7}{e, x) 

where x is a fixed, smooth cutoff function with support on {(x, y) : y < rj(e, x) + 2} and 
chosen so that T(e, •) is a diffeomorphism. Because we are fixing a representation for 
T, the only unknowns in the problem are ipj and Q- In particular, T is determined 
entirely by 77. We remark that this map has been the basis for a number of studies of 
water waves. To reference only the most immediately relevant results, we point out 
that Wahlen, Ehrstrom et al. rely on it in their investigations of steady waves with 
critical layers (cf. [21, 11, 10]). 

For notational convenience, we denote the inverse of T(e, •) by S(e, •). We shall 
also use the convention that the coordinates in Q are in the variables (x,y), while 
the unbarred variables indicates coordinates in f2(e), e =/= 0. 

Define the transformed stream function "J/ by the relation 

1>(e r ) = [9(e,-)]oT(e,-) ) 

or, equivalently, 

*(£,■) ~Me,-)}oS(e,-). (5.5) 

We shall suppose that *&(Q,x,y) — ^ (y), meaning that the unperturbed flow is 
laminar. 

Let di denote partial differentiation with respect to the i-th physical variable, for 
i = 1,2, and, for any function / of two variables, let f ± := dif. Then an elementary 
computation confirms that, —Atp(e,-) — in fi(e) if and only if ^(e, •) satisfies 

the following equation in Qq: 

£{r))^ := A l0 did 3 ^ + B t d^ + 7(f) = in fi , for all e. (5.0) 

Here we are adopting the summation convention over repeated indices, 7 represents 
the vorticity strength function for the flow (7 = for irrotational flow), and 

An = Av{n) := [{dkXKe, •)(^T J -)(e, •)] ° S(e, •) 

Note that are we stating that £ depends on 77, while only T occurs above. This is 
valid because T is determined uniquely by 77 in view of (5.4). 

The jump condition on the boundary in (5.1) can likewise be reformulated in 
terms of ^, resulting in the following 

[(C^*) 2 ] + 2g [p] H - Q = on y = 0, (5.8) 

where H(e,x,y) := S2(e,x,y) — ^(e, x, — 1) is the height above the ocean bed in the 
terms of the coordinates (x, y), and 

C ii =C«(T 7 ):=[(a j -T i )(e J -)]o5(e,-). (5.9) 
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The Dirichlet conditions for ^ derive from those for ip and the definition of T: 

! = n P ° ° ng - = ~ h (5.10) 

* = on y = 0. K J 

Lastly, the Neumann boundary condition at y = +oo translates to the same 
condition for since we have that T asymptotically approaches the identity. 

V" L, J / — > (—A, 0) as y — > oo uniformly in x and e. (5-H) 

In summary, we find that there exists a non-laminar solution to (ip,r],Q) to 
(5.1) provided that there is exists (r],X,Q) for which there are nontrivial solutions 
('J, H, Q(A)) to (5.6)-(5.11). To make the latter problem tractable, we suppose that 
there is a unique solution \E' for given (77, A). This is valid since, when 77 = 0, T = lq 
and thus £(A,0) = A. It follows that the operator is an isomorphism for e in a 
neighborhood of 0. Explicitly, we define — ^^'(A, 77) to be the solution of 

£(A,T7)#W=0 in^ 1} , 

*W =0 on y = 0, (5.12) 

y-'-'I'W — > (— A,0) as y — > 00 uniformly in x. 

and define 4 , ^ 2 - ) = ^^'(A, rf) to be the unique solution of 

£(A,T7)tf( 2 > =0 in^ 2) , 

*W =0 on y = 0, (5.13) 

\|/( 2 ) = — p Q on y — _1. 

Then solving (5.6)— (5. 1 1) is equivalent to the following: find (A, 77) such that G(\, rj) = 
0, where 

g(\,rj) := [(ffyfa.A^fa.A)) 2 ] +2 ff W( J ff(t ? ))| fi= o-0. (5-14) 

Notice that a laminar flow corresponds to a solution where 77 = 0, and if = 7r 2 + 1 
(since we have dictated that T\ = tti, this is equivalent to saying T = l). 

5.2. Linearization. Echoing the approach in §3 and §4, we begin by proving 
the existence of a one-parameter family of laminar flows. We the proceed to linearize 
(5.14) along this family, in order to lay the groundwork for a bifurcation theory 
argument. 

Lemma 5.2 (Laminar flows). There exists a one-parameter family of laminar so- 
lutions (Q(A), 0) to (5.14) for A > 0. The corresponding family of laminar transformed 
streamed function ^o(X) are given by 

*o(A, y ) = ( ~ X ? f f ° r y>°> (5.15) 
K yl \ Pay for y < 0, y ' 

and 



Q(A) := 2(A 2 -pg) + 2. 9 M 



(5.16) 
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Proof. For the flow to be laminar, we must have that H = ~k-i + 1, so that, from 
(5.1) we see that the transformed laminar stream function vPq = ^ / o(^j 7r 2) is harmonic 
in SIq 1 - 1 U fig 2 ' and has the following boundary data: 

#o(0)=0, # (-l) = -po, *o -> ~ x ^ V -> °°. 

We deduce that it must be as in (5.15). The formula for Q(X) then follows from 
(5.14). □ 

Fixing A, we now compute the Frechet derivatives of E and Q . In what follows, D v 
denotes Frechet differentiation with respect to 77, while d e is the (finite-dimensional) 
partial derivative with respect to e. Also, where there is no risk of confusion, we shall 
suppress the A dependence. 

Let the variations be denoted by 

h := (d e H)U=o, $ := (W| e =o, <\> ■= (^)| £ =o, 
C := (d tV )\ £=0 , r:=(9 e T)| e=0 , a := (d e S)\ e=0 . 

Remark 6. There are a few points that should be made here: 

(i) (f> is not continuous over the x-axis (indeed, it is not well-defined there at 

all). 

(ii) An elementary calculus results states r = —ex. 

(iii) Chasing the definitions, it is obvious that a% = h. 

(iv) Since \P vanishes on the x-axis for all values of e, we have 

= ^(e,x,T 2 (e,x,T][e,x))) 
= <f(x, 0) + *' (0) (r 2 + T 2 , 2 (0, x, 0)C(a?)) 

By (ii) and the fact that T 2 (0, •) = tt 2 (-)j this implies 

h = C on V — 0. 

in other words, £ is t/ie irace 0//1 on i/ie x-axis. 
Now 

eJ 4y = [T,, fc£ (e, -)^, fc (e, •) + T iifc (e, •)^, fce (e, ■)] ^fo 

+ S e . e (e, -)[T iM (e, -)Tj tk (e, •) + T jM (e, -)T t , k {e, •)] o 5(e, •), 

and thus, evaluating at e = we find 

(AjAjjXO) = o"j fe Ti,fe + ^Tj,fc = diTj + d 3 Ti. 

Here we have used the fact that 5(0, •) = T(0, •) = lq . 

The calculation of the linearization of the first-order coefficients proceeds in the 
same fashion: 

d e B t = [T idje (e, •)] o 5(e, •) + 5 M (e, -)[T iJjk (e, •)] o 5(e, •) 
D^ i (0) = r i ,^+o- fc [^ ifc ] 
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Collecting these two facts, we see that 

(£„(A, 0), *) = Ay (0)0*9,$ + S<(0)flSi* 

= A$-%Ah. (5.17) 

Note that the last line follows from observing that t 2 = —cr 2 = —h. 

Remark 7. This formula can also be obtained formally by noting that 

d e Me,-)]=d e Me,S{e,-))} 

= te(e, S(e, •)) + ^x(e, S(e, ■))5 1 , e (e, •) + ^ 2 (e, S(e, ■))S 2 ^e, •), 

and £/ius, evaluating at e = 0, we have 

$ = 4> + %h. 

Taking the Laplacian of both sides of the equation leads to (5.17). 

Next we consider the linearization of Q, the transmission boundary condition. 

d e C l3 = [T itje (e, •)] o S(e, •) + S fe , e (e, ■)[T iJk (e, •)] o S(e, •), 

whence 

(D v Cij)(0) = T id + a k d k (Sij) = -dj(Ji. 
Using this identity with (5.14), we compute 

(0„(O,A),C) =2l({(D ri C ij )(0),h)d i ^ Q + C ij (0)d^)(C ij (0)d i ^ )j 

+ 2sMC 

= 2 l(-%d 2 h + d2*)*' j + 2g [p] C- (5.18) 
As a consequence of our choice of T in (5.4), we have 

S 2 (e, x, y) = (?7(e, x) + l)y + rj(e, x) in a neighborhood of {y = 0}, 
and thus 

h(x, y) = C,(x)y + i n a neighborhood of {y = 0}. 

This implies = C 011 2/ — 0- Therefore (5.18) can be written 

(a,(A, o), o = 2 ( 5 m - 1(*^) 2 ] ) c + 2 i%a 2 n . 

Finally, to make sense of this, we compute explicitly the dependence of d 2 § on £. 
From (5.17) we have that $ solves 

A$ = %Ah in Q \ {y = 0}, 

° n ^ = ' 1 (5-19) 
$ = onj/= — 1, v ; 

V-L$(!) -> (0,0) as y —¥ oo uniformly in x and e. 
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Equivalently: 

f A(¥{,* - (%) 2 h) = mft \{y = 0}, 

(%$ - (%) 2 h)W = -C((*' ) 2 ) W on y = 0, 

(^$-(*[,) 2 /i)( 2 ) =-C((*o) 2 ) (2) ony = 0, (5.20) 

%<!>-(%) 2 h = ony = -l, 

V J -(^E'o < J' _ (^o) 2 ' 1 ) — * (0' 0) as y — ► oo uniformly in x and e. 

Taking T := — ($>' ) 2 h, (5.20) makes clear that T depends only on the trace of 
h on {y = 0}, i.e., on £. Moreover, because of the evenness and periodicity, we can 
compute T = T(£) quite explicitly. Letting • designate the Fourier transform in the 
x-coordinate with Fourier variable k, we have 



(5.21) 



(5.23) 



(TW)(fc,y) = -A 2 C(fc) (cosh(fcy) - sinh(%)) 
(T^){k, y) = -poC(fc) (coth (k) sinh (ky) + cosh (ky)) 
From this we readily obtain 

[&?!(&) = K(k) (p 2 cothk- A 2 ) = (p 2 L»cothD- A 2 L>)C. (5.22) 
That is, we may view Gr){X,0) as a Fourier multiplier: 

(6,(A,0),0=2 S WC + 2W] 

= 2(p 2 D coth D - X 2 D + g 

Here (abusing notation slightly) D = —id s - 

5.3. Proof of local bifurcation. With the expressions for Q n (X,0) derived in 
the previous subsection, the proof of the local bifurcation theorem is relatively simple. 
Define 

X ■= C 2 p + a ([-7T,7r}), Y := C*+°<([-7r,7r]). 

In what follows, Q is considered as an operator with domain M x X and codomain Y. 

Lemma 5.3 (Null space and range). Under the assumption that the local bifur- 
cation condition (5.2) holds, there exists a A* > such that the following statements 
are true. 

(i) Q v ([\*,0) is a Fredholm operator of index 0; 

(ii) A/"(^(A*,0)) is one- dimensional and spanned by C (x) := cosai; moreover 

(iii) £ € n(Gr,(X*,0)) if and only = 0. 

Proof. By (5.23), we know that C € M(G V (X* , 0)) if and only if it satisfies 
(p 2 fccothfc - X 2 k + g M)C(fc) = 0, for all k > 0. 

For each A > 0, 

m(k; A) := p\k coth k — X 2 k + g fp} 

is injective as a function of k, and, fixing k > 0, m(fe; A) — > — oo asA-> oo. The local 
bifurcation condition implies that m(l; A) > 0, and hence there exists a A* such that 

plkcothk- (\*) 2 k + gfpj =0 if and only if k = 1. (5.24) 
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For this choice of A we must have ((k) = for k ^ 1, and thus evenness dictates that 
C G span(cos (x)). We conclude the null space is one-dimensional, proving (b). 
On the other hand, suppose that £ € lZ(G n (A*, 0)). Then 

i(k) = m(k;X*)C(k), for k > 0, 

and, in particular, £(1) = 0, by the definition of A* in (5.24). It follows that £(1) = 
is a necessary condition for inclusion in the range. Conversely, if £ is any element of 
Y with |(1) = 0, then (el defined by 

is in the preimage of £ under G T] (X* , 0). Note that this definition is permissible since 
m(k;X*) is nonvanishing for k ^ 1, again by (5.24). We have therefore shown that 
6 € 7Z(G V (X*, 0)) if and only if £(1) = 0, which is (iii). Of course, this implies 
immediately that the codimension of the range is one, and so (i) follows. □ 

Remark 8. Equation (5.24) in fact gives an explicit definition for the point of 
bifurcation: 



X* = J p 2 cothl + g {pi- 



p-roof of Theorem 5.1. In the previous lemma, we confirmed hypothesis (iii) of 
Theorem A.l; while (i) and (ii) clearly hold. All that remains is the transversality 
condition, hypothesis (iv). But observe that 

({g n x(X*,0),C)f(k) = -4A*fcC(fc) = -4A*£ U . 

By Lemma 5.23 (c), (Gr)\(X*, 0), £*) cannot be an element of the range of G v (X*,0). 
The statement of the theorem then follows from a straightforward application of 
Theorem A.l. □ 

6. Local bifurcation for shear flow in the atmosphere. In fact, a fairly 
simple extension of Theorem 5.1 is possible when we assume that the flow in the 
atmosphere region has constant vorticity. Let 

_ / 70 in f>W 
7 - \ in 

where 70 is a fixed constant. Rather than study the relative stream function directly, 
we instead look at the perturbation of the stream function for the background shear 
flow. In other words, let 

so that 

= Atp + 7 = 0, in fl. 

We shall call ip the modified stream function. The boundary conditions for ip follow 
naturally from those for tp enumerated in (5.1). Notice, however, that the interpre- 
tation of the condition V ip — > (~X, 0) as y — > 00 is different: we are now imposing 
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the precise way in which the shear of the velocity field approaches infinity, instead of 
the limiting value of the velocity. Even more importantly, we point out that ip is not 
continuous due to the jump in 7 over the interface. The appropriate choice of Banach 
spaces in this setting is therefore 

y : = v , Qi v) € y '■ < c, (u — Uy,v) —> (0, 0) uniformly as y —> 00} 
:= {(Q, ^ r}) : (Q, V, v) € S" , -> (-A, 0) uniformly as y -> 00}. 

With the notation established, we can now state our main result for the shear 
flow. 

Theorem 6.1 (Local bifurcation for unbounded shear atmosphere). Suppose the 
following local bifurcation condition holds 

^cothl + «?M-i(7 -) 2 + i 7 o7o~ >0, (6.1) 

where 7^ := min{7o,0}, then the following statements hold. 

(i) There exists a continuous curve of non-laminar solutions to the stream func- 
tion equation for irrotational flow in the air and water, and unbounded atmosphere 
region (5.1) 

C[ oc = {(Q(s)Jj( S ),n( S ))e^:\s\<e}, 

for e > sufficiently small, such that (Q(0),ip(0)) = (Q*, 1 ^^), and, in a sufficiently 
small neighborhood of (Q*,^) in K x X, C[ oc comprises all non-laminar solutions. 
In particular, 

v (s)(-) = s cos (-)+o(s) inC%£*(R). 

(ii) There is a corresponding continuous curve 

Cloc - {(Q(s), u(s), v(s), Q{s),P{s), 17(a)) € ,^0 ■ \s\ < e} 

of small amplitude solutions to the Eulerian problem which likewise captures all non- 
laminar solutions in a sufficiently small neighborhood of the point of bifurcation. 

Remark 9. As in Theorem 5.1, a simple generalization of this result is possible 
we have depth d > and allow perturbations of the laminar flow with minimal period 
2irk* , where k* is the smallest nonnegative integer satisfying 

^k* coth (k*d) + g Ipj - i( 7o -) 2 + i 7 o7o~ > 0- 

We also note that when 70 > 0, (6.1) reduces to (5.2). 

6.1. The transformed problem. Define the transformations T as in §5.1, and 

put 

$ := ifo S. 

Since ip is harmonic in the unknown domain, \& will be in the kernel of the elliptic 
operator £(rf), just as $ was in the previous section. The main difference will come 
in the boundary condition, and, most significantly, in the transmission boundary 
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condition. Writing ^ — 7y 2 /2 = \& and inserting this in to the Bernoulli equation, we 
find that the entire problem is equivalent to the vanishing of the operator 



G(X, V ):= (C i3 (v)d^( V , A)) 1 



+ 2sW (H( V ))\y =0 -Q, 
where ^(A,^) is defined to be the (unique) solution to 



(6.2) 



E(i})V = 

~ y s l =° on y = 



in n \{y = 0} 

on y = 
on y — 



$(2) = 

$ = — po on y = — 1 

V- L $^(-A,0) asy^oo. 



(6.3) 



The laminar solution = ^(A, 0) is the same as for the ideal flow, i.e. \l/o = ^o- 
This is simply because the only way in which the vorticity appears in (6.3) is as a 
coefficient of 52, but 52(0) = 0. We record this fact in the following lemma. 

Lemma 6.2 (Laminar flows). There exists a one-parameter family of laminar 
solutions (Q(A),0) to (5.14) for A > 0. The corresponding family of modified laminar 
transformed streamed function ^o(X) are given by 



*o(A,y) = 



-Xy for y > 0, 
PoV for y < 0, 



(6.4) 



Q(A) :=2(X 2 -p 2 ) + 2g\ 



(6.5) 



Remark 10. Recalling the definition of the modified stream function, we have 

*o = *o + |y 2 , 
where ^o is the stream function for the laminar flow. Thus, 

% = joy = -X- joy- 
In other words, if 70 < 0, then there is a critical layer in the air at y = |7o|/A. 
Of course, because \&g" vanishes identically, this will be a neutered layer in the sense 
discussed in the introduction. 

6.2. Proof of local bifurcation. Next consider the linearization of the opera- 
tors £ and Q around the laminar flows. Let the variation of \& be denoted by <&, and 
otherwise adopt the same notation as in the ideal atmosphere case considered in§5. 
Then 



<S,(A,0),C> 



(-^ d 2 h + d 2 <Z>)% - 2 7o CA + 2.g M C 



2 UM- 



(%) 2 -2 7o A C + 2 %d 2 <5> 



(6.6) 
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On the other hand, solves (5.19), and so we may define T := ty' $ — (^q) 2 /i, and 
from an identical argument we find 



d 2 T 



which can be understood as a Fourier multiplier using (5.22): 
(<7„(A,0),0 = 2(p 2 DcotlaD-X 2 D+gl 



7oA)C- 



(6.7) 



Observe that the only difference between the symbol in (6.7) and that in (5.23) is the 
— 27oA£ term. 

Lemma 6.3 (Null space and range). Assume that the local bifurcation condition 
holds. Then there exists a A* > such that the following statements hold: 

(i) S v (X*,0) is a Fredholm operator of index 0; 

(ii) M{Qn{X* , 0)) is one- dimensional and spanned by C*(x) := cosx; moreover 
(hi) £ € n(G n (X*,0)) if and only = 0. 

Proof. From (6.7) it is clear that £ £ M(Q V (X* , 0)) if and only if 



m(k;\*)((k) = 0, 



for all k > 0, 



where 



m(k; A) := pgfccothfc - A 2 fc + g [p] - 7 A. 

First we note that, for any fixed A > 0, fc i— > m(fc; A) is a strictly decreasing function 
and thus any root is unique. Since we are primarily interested in the case where the 
perturbations of the laminar flow are minimally 27r-periodic, we wish to consider A 
for which m(l;A) = 0. Differentiating, it becomes clear that m(l; A) is a concave 
function of A tending to — oo as A — > oo and that the maximum occurs at 



A = (-?) 



1 



;7 • 



The local bifurcation condition (6.1) implies that 

maxm(l; A) > m(l; — -7rT) > 0. 

By continuity, we have that there exists a A* such that m(l; A*) = and fh(k; A*) 7^ 
for k 7^ 1. This proves (ii). The proof of the remaining parts follows exactly as in 
Lemma 5.3. □ 

Proof of Theorem 6.1. We have already laid the groundwork for a bifurcation 
argument via Crandall-Rabinowitz. The only detail that remains to be checked is the 
transversality condition. With that in mind, we calculate 

((0,a(A*,0),O)» = -4A*fc<f*(*0 ~ 2 70 <f*(fc) 
= -4A*<5ifc - 2j Si k - 



In light of Lemma 6.3 (c), the only way for C* to be an element of 1Z(G 71 (0, A*)) is if 
2A* = —70- But this scenario is excluded by the local bifurcation condition (6.1), as 
m(l; —7^/2) > while m(l; A*) = 0. This confirms that the transversality condition 
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holds, and thus we obtain the theorem via a routine application of Theorem A.l. 
□ 

Appendix. Here we present two standard theorems. The first is the classical 
work of Crandall-Rabinowitz on bifurcation from simple (generalized) eigenvalues: 

Theorem A.l. (Crandall-Rabinowitz [8]) Let X andY be Banach spaces, I C K 
an open interval with A* G /. Suppose that T : I x X — » Y is a continuous map with 
the following properties: 

(i) T(\, 0) = 0, for all A G I; 

(ii) D\T ' , D2J- and D1D2J 7 exist and are continuous, where Di denotes the 
Frechet derivative with respect to the i-th coordinate; 

(iii) D2J-(\*,0) is a Fredholm operator of index 0, in particular, the null space 
is one- dimensional and spanned by some element w*; 

(iv) D 1 D 2 F{K,Q)w,i'R{D 2 F{K 1 Q)). 

Then there exists a continuous local bifurcation curve {(A(s), w(s)) G M x X : \s\ < e} 
with e > sufficiently small such that (A(0), w(0)) — (\#,w t ), and 

{(A,w) eU 0,T(X,w) = 0} = {(X(s),w(s)) e K x X : \s\ < e} 

for some neighborhood hi of (A*, 0) in R X X . Moreover, we have 

w(s) = sw* + o(s), in X, |s| < e. 

If D\T exists and is continuous, then the curve is of class C . The reader should 
view this as a special case of the more general Lyapunov-Schmidt reduction procedure, 
a good discussion of which can be found in [6] . 

The second theorem is on the Fredholm solvability and regularity of solutions 
to linear elliptic equations with jump conditions across an interface. A classical ref- 
erence for this is [14]. We quote here a parsed version of Theorem 16.1 from that 
book, incorporating the discussion preceding and following the theorem statement 
and simplifying the hypotheses to better match our needs. 

Let = £1^ U fl^ be a domain in K 2 with smooth boundary of class C 2+a , 
where and 51^ 2 ^ are the connected components of fl and their shared boundary 
I := dfl^ n <9f^ 2 ) is a simple curve (which will also be of class C 2+a ). Consider the 
elliptic problem 

d Xi iflijd Xj u) + bid Xi u + cu = f in Q 

H = on 1 r AA ) 

ldd u u} + au = g on I 

u = h on <9fi \ I. 

Here we are using summation conventions, d v denotes the conormal derivative with 
respect to a^, a is a real constant, and the coefficients and forcing terms are assumed 
to have the following regularity: 

aij uniformly elliptic, d > do > 0, 

h, c, f g c a (nW)nc a (oW), 

a ih d, gec 1+a (nw)nc 1+a (n( 2 )), [ ' 

fteC 2+a (J]W)nC 2+a (!]( 2 )). 

Theorem A. 2. (Ladyzhenskaya-Ural'tseva [14, Theorem 16.1]) Consider the 
elliptic problem (A.l) with the assumptions listed in (A. 2). Then the following state- 
ments hold. 
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(i) (Frcdholm solvability) The existence of a classical solution 



u g c 2+a (n \i) n c 1+a (n<U) n c 1+a (nW) 

is implied by the uniqueness of solutions. 

(ii) (Elliptic regularity) Note that, due to the fact that the [u] = 0, a classical 
solution can be extended to £1 as a continuous function. In fact, if a classical solution 
u exists, then its extension to f2 is of class C a (Q). 
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